Abstract 



Double, Glen Paul. Investigations of Lepton and Baryon Acceleration in Relativistic 
Astrophysical Shocks. (Under the direction of Donald C. Ellison.) 

Gamma-ray bursts, occurring randomly each day anywhere in the universe, may 
be the brightest objects in the sky during their short life. Particle acceleration in 
mildly relativistic shocks, internal to the main blastwave, may explain the early in- 
tensity peaks in gamma-ray bursts and the afterglow may be explained by energetic 
particles accelerated by the main ultrarelativistic blastwave shock as it slows to the 
mildly relativistic range. To help explain the phenomena, a nonlinear relativistic 
Monte Carlo model was developed and used to study lepton and baryon acceleration 
by mildly relativistic modified shocks with the magnetic field parallel to the shock 
normal. The study showed that for equal densities of leptons and baryons, lepton 
acceleration is highly sensitive to the shock velocity profile. With the shock fully 
modified by energetic baryons, the injection efficiency of leptons, relative to baryons, 
increases with Lorentz factor, and injection efficiency will reach a maximum well be- 
low that of baryons at the same momentum. Given the assumptions in this model, 
if the particles are energized by shock acceleration, leptons will always carry far less 
energy than baryons when the lepton and baryon densities are of the same order. 
It was determined that the lepton to baryon number density ratio must be approx- 
imately 3 X 10^ for both species to equally share the kinetic energy of the shock. 
This energy equipartition density ratio is independent of shock speed over the range 
of mildly relativistic Lorentz factors used in the study, but the result may extend 
to ultrarelativistic speeds. The study was a special case of a larger effort that will 
include relativistic oblique modified shocks and computer generated gamma-ray spec- 
tra when the model is completed. New magnetohydrodynamic conservation laws and 
relativistic jump conditions were developed for the model, along with a new equation 
of state and a new method for estimating the adiabatic index in the mildly relativistic 
range. The present state of the model shows smooth transitions of shock parameters 
from nonrelativistic to highly relativistic unmodified shocks while allowing oblique 
magnetic fields and a pressure tensor. 
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Chapter 1 



Introduction 



Shock waves occur in fluids whenever a disturbance propagates faster than the characteristic speed 
of sound in the fluid. The same phenomena occur in space where the fluid is a plasma consisting 
of positively and negatively charged particles permeated by magnetic fields and negligible electric 
fields. 

In space where the particle density is low, typically a few particles per cubic centimeter, shock 
waves are collisionless rather than acoustic-ion shocks because the effective diffusion length due to 
charged particles interacting with the magnetic field is generally much shorter than the mean free 
path between particles. An acoustic-ion shock wave propagates by Coulomb forces, for example 
shock waves in the Earth's atmosphere, but collisionless shocks propagate through disturbances in 
the magnetic field of the space plasma. 

Shock waves are found throughout the universe. Examples of nonrelativistic shocks range from 
transient disturbances on the surface of the sun ( HoUwcg , 1982 ), a standing bowshock from the solar 
wind blowing across t he Earth's magnetic field tellisoJ . Il985() . to energetic shocks associated with 
supernova explosions i Dwarkadas fc ChevalieIlll998^ and other extreme types of exploding objects. 

Relativistic shocks originate from much more energetic events, including the most powerful 
supernova explosions fhvDernovae: lR.eesl . l200f]l) . jets in extragalactic radio sources, and the mysterious 
objects that produce gamma-ray bursts. The study of astrophysical shock waves is important because 
the shock waves can accelerate charged particles, by first-order Fermi acceleration, to very high 
energies and possibly explain the existence of the highest energy cosmic rays. The study of relativistic 
shocks may bring some understanding to the spectra seen in gamma-ray bursts. 

As an introduction, an overview will be presented of gamma-ray bursts, diffusive shock accelera- 
tion and the mathematical and computational tools that were developed by this author and others to 
study particle acceleration by relativistic shocks. The introduction is concluded with the objectives 
and the overall organization of this dissertation. 



1.1 Gamma-ray Bursts 

Gamma-ray bursts were discovered accidentally in the 1960's by satellites intended for monitoring 
nuclear testing. Since then there has been an abundance of research, both observational and the- 
oretical, in an attempt to better understand the nature and origin of the bursts. Much has been 
published and an excellent review, among a number of others, was published by lRccs ( 20q3)- 

Gamma-ray bursts are extremely energetic blasts (~ 10^^ ergs if the blast is isotropic; ~ 10^^ 
ergs if beamed) of primarily gamma and X-ray radiation that are often the brightest objects in the 
sky. They occur daily and seem to come from every direction, from our galaxy to near the edge of 
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the observable universe. Gamma-ray bursts remain a great mystery, and only recently were optical 
counterparts to the gamma-ray bursts observed. Virtually nothing is known about the engine that 
drives the burst. An understanding of gamma-ray bursts will lead to a better understanding of the 
form and evolution of the universe as a whole. 

1.1.1 Observational Overview 

Gamma-ray bursts typically have a duration of seconds to minutes, as shown in Figure lTTTl but some 
have afterglows lasting weeks or months. The bursts are isotropic on the sky and, for those which 
have been associated with host galaxies, show large red shifts which imply cosmological distances, 
and therefore extremely large energy output to produce the brightness we see. Most of the power 
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Figure 1.1: Typical time plot of a gamma ray burst. Longer bursts tend to have a more complex 
structure. 

of the burst is non-thermal radiation emitted in the 100 - 1000 keV range, usually with a number 
of short duration intensity spikes which imply a compact size with a radius not more than a few 
thousand kilometers. A "typical" energy spectrum is shown in Figure ^21 i-e., typical in the sense 
that The burst normally has a brief period of intense multiple energy peaks, but every one is different 
and it makes it difficult to summarize their basic features. The radiation energy is usually above 
50 keV, sometimes in the MeV range, but many gamma-ray bursts have also exhibited an optical 
afterglow that follows the initial burst, some with energies dipping into the infrared and microwave 
regions. The expansion rates of gamma-ray bursts are explained by ultrarelativistic flows with 
Lorentz factors greater than 100, possibly up to 1000 or more. 

1.1.2 Proposed Models 

The actual inner engine that produces the gamma-ray burst is unknown. Researchers speculate that 
compact objects (black holes, neutron stars, etc.) somehow merge and produce the ene rgy to power 
a gamma-ray burst. The most widely accepted model appears to be the Fireball model l|Piranl 11999') 
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Figure 1.2: Typical energy spectral plot of a gamma-ray burst. 



shown in Figm'e [T31 By some unknown means, a gamma-ray and/or neutrino radiation-dominated 
highly relativistic blastwave gives rise to an external shock and an electron-positron pair-dominated 
ultrarelativistic wind in a region of low baryon density. A reverse shock, traveling backward into the 
lepton wind, and possibly internal shocks interacting with each other or the reverse shock at slower 
but still relativistic speeds, may be responsible for the series of intensity peaks seen in the initial 
gamma-ray burst. The forward shock slows quickly (on the order of minutes) to mildly relativistic 
speeds as it picks up ambient mass in the interstellar medium. Kinetic energy in the outflow is 
randomized by magnetic fields, probably collisionless shocks, and is converted back into radiation 
seen as the afterglow. The computational methods used here allow the study of lepton and baryon 
dynamics (i.e., Fermi acceleration and the ensuing nonthermal distributions of energetic particles) in 
the internal shock regions and the external afterglow regions where the shocks speeds have Lorentz 
factors of 20 or less. 



1.2 First order Fermi acceleration 

Acceleration of a c harged particle as a result of diffusive scattering in converging flows was first 
suggested bv iFerrni (1949). There were a number of important papers that developed the con- 
cept of diffusive acceleration, or the steady diffusive acceleration of test particles as they cross the 
shock boundary numerous times and scatter through stochastic processes. Due to the random na- 
ture of the scattering process, information is destroyed (i.e., entropy is increased) and the resulting 
energetic particle distributio n is relatively i ndependent of t h e details of the particle energy d istri- 
bution prior to acceleration l|Drurvl Il983|) . lKrvmsk\l l)l977(l . IXxford. Leer, fc SkadronI l)l977j) . and 
[Blandford & Ostriker ( 1978) used a macroscopic approach to develop an analytical theory of particle 
acceleration, assuming isotropy in the particle momentum distribution in the fluid frame (because 
the particle velocity is assumed to be much higher than the fluid velocity), and assuming collision- 
less, elastic scattering, and assuming the scattering centers are frozen into the background plasma. 
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forward shock 




Figure 1.3: The Fireball Model. A highly relativistic (quite likely beamed) forward shock moves 
outward into the ISM and a reverse shock moves back into the outward moving lepton wind. Various 
relativistic internal shocks may interact with each other or with the reverse shock to produce intense 
radiation spikes that make up the initial gamma-ray burst. 



An accelerated particle momentum spectrum can be generated by using the transport equation in 
the local fluid frame, and matching boundary conditions for the steady-state solution on each side of 
the infinite plane shock. This approach is straightforward, but it does not provide an understanding 
of the physical processes that are responsible for th e par t icle acc e leration. A micr oscopic de rivation 
of diffusive particle acceleration was presented by iBelll ( 1978a^ , IPeacockl l(l98l|) , and also iMichell 



which provi des much mo r e insight into the physics of the p article acceleration pro cess. A 
number of papers, iDrurvl ||1983|) . iBlandford fc Eichleil ()l987l) . and I Jones fc EllisonI l)l99l|) among 
others, review both approaches in detail, covering the effects of oblique magnetic fields, modification 
of the the shock velocity profile by the backpressure of accelerated particles, and the connection 
to cosmic rays. The papers cited above deal with nonrelativistic shocks, but they provide a good 
introduction to shock acceleration of particles. 

The most basic concept of nonrelativistic first order Fermi acceleration is shown below. A strong 
shock wave provides the mechanism by which particles are scattered in pitch angle by magnetic wave 
turbulence on each side of the shock, as shown in Figure [TM Three frames of reference are shown. 
The upstream plasma frame sees the shock frame moving toward it with velocity —uq and the 
downstream frame moving toward it with relative velocity Urci = U2 ~ uq. The shock frame sees 
the upstream frame as a flow moving toward it with a velocity -|-mo and the downstream moving 
away with a velocity +U2 The downstream frame sees the upstream frame moving toward it with 
relative velocity Wrci — uq~u2 (i.e., the upstream and downstream frames are converging at a speed 
Mrei) and the downstream frame sees the shock moving away and to the left with a velocity — U2- 
A charged particle, for example a proton, initially at rest in the upstream frame crosses the shock 
from upstream to downstream because the shock is moving toward it. In the downstream frame, the 
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Figure 1.4: Basic conceptional diagram of a plane shock moving through space with a particle 
scattering elastically in the downstream frame and gaining energy. For nonrelativistic shocks, the 
upstream and downstream frames converge with relative velocity Urci = uq — U2- 
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particle has velocity +Mic1- Suppose the particle elastically scatters off of magnetic turbulence in the 
downstream frame as shown in Figure IT"^ and ends up moving to the left with velocity — Urci- As 
the particle moves back across the shock, the upstream frame sees the particle, which was originally 
at rest, now with a velocity of — 2mioi (because the particle had a final velocity of —Mid in the 
downstream frame and the downstream frame is moving toward the upstream frame with a velocity 
of — Uiei- Hence, the particle gained in kinetic energy by making a complete passage across the shock 
and back into the upstream frame. The same argument applies to a particle originally downstream; 
the particle will also gain the same amount of kinetic energy. By making numerous passages back 
and forth across the shock, a particle can gain a large amount of kinetic energy, easily becoming 
relativistic. Obviously, most particles are at some angle with respect to the shock normal, therefore 
only the x component of the particle's velocity contributes to the particle's gain in energy. Since the 
process is completely random, a distribution of particle ene rgies will result that is independent of 
the original state of the upstream conditions r ilDrurvllTQsi . In fact, the distribution will result in 
a power law whose slope depends only on the compression ratio, r = uq/u2, between the unshocked, 
upstream plasma and the shocked downstream plasma, the details of which are discussed at length 
in the papers cited above. 

The main points of the derivation of the power law for accelerated particles are described in 
a review bv IPrury (1983). Suppose a particle has momentum p, velocity v and pitch fi = cos{(f)) 
between the momentum vector and the x axis; then the average flux weighted change in a particle's 
momentum with respect to the local plasma frame when it crosses the shock is 



(Ap) = p I [^i{uo - U2)/w]2/id^ gP(wo - U2)/v 



(1.1) 



Therefore, pi ~ po[l + (Ap)], and after N shock crossings (or N/2 returns from downstream), the 
average momentum will be 



w 2 4 

(Pn) = + ^(mo - U2)/v,]po W ]^[1 + -{Uo - U2)/Vi\pQ 



(1.2) 



Then 



In 



PN 



PO 



4 "^^^ 1 

-{U0-U2)Y.- 



(1.3) 



From the downstream frame, the flux of particles passing through the shock from downstream to 
upstream is 



n {u2+ v^)dv^ 



\{u2~vf 



and the flux of particles passing through the shock from upstream to downstream is 



(1.4) 



n / {u2 + Vx)dvx^ -{u2+vY 



(1.5) 



The probability of return from downstream to upstream is the ratio of the two fluxes, or 



Pr = 



/ U2-V 
\U2 +V 



(1.6) 
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The probability that a particle has returned to the shock at least N/2 times is 



Pr{N) 



N/2 ^ 
' U2- Vi 



n 



U2 + Vi 



(1.7) 



Next, equate summations from the probability equation above and the previous momentum equation: 

\n[Pnm^-^^^n(P-) (1.8) 



Pnip) -(f) (1-9) 



This gives the probability that a particle will reach at least momentum p: 

-3m2/(mo— M2) 

.Pa 

The number density of particles accelerated to momentum p is the product of the initial number 
density and the probability function: 

nip)=Pn{N)n{po)^no^{^) (1.10) 
Finally, the distribution function is 

Therefore, with the assumptions of test particles (which leave the plane shock structure unmodified), 
particle velocities far greater than the nonrelativistic shock velocity and an isotropic distribution of 
particle momenta, the resulting momentum distribution spectrum depends only on the compression 
ratio r, i.e., 

f{p)ocp~"' where ^ r~~T (l-^^) 

Although the derivation of the functional relationship for the momentum distribution with its 
accompanying assumptions were focused on nonrelativistic shocks, the same concepts of particle 
acceleration can be applied to relativistic shocks, provided that the new properties which apply to 
relativistic shocks are taken into account. For example, particle velocities are close to those of the 
shock velocity (i.e., the speed of light), and this leads to different characteristics and criteria for the 
particles interacting with the shock. The probability of return equation H1.6|l is true in general, even 
at relativistic speeds, provided the downstream momentum distribution is isotropic, and it should 
always be found to be isotropic at least one or two mean free paths downstream from the shock. 
However, the momentum distribution at the shock is highly anisotropic. Particle acceleration by 
relativistic shocks will be addressed in a later chapter. 



1.3 Conservation laws and jump conditions 

Particles must obey conservation of momentum, energy and particle continuity across the shock, as 
well as obey Maxwell's equations and the thermod ynamic equation of state. T he nonrelativistic jump 
conditions are well known; for example, in .Ellison. Baring, fc JonesI l|l996l) . the jump conditions 
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for nonrelativistic shocks with magnetic fields at oblique angles are given in full, and implicitly 
incorporate the adiabatic equation of state. The jump conditions relate the upstream unshocked 
plasma to the downstream shocked plasma. Hence, for particles scattering back and forth across the 
shock, the jump conditions must be satisfied at each crossing and provide a self-consistent solution. 
This topic will be generalized for relativistic shocks and will be thoroughly discussed in the next two 
chapters. 

1.4 The Monte Carlo approach 

Particle acceleration by shocks is inherently complicated and leads to nonlinear differential equations. 
For example, the compression ratio depends on the adiabatic index, wh i ch is affected by the same 
shock speed that determines the compression ratio teUis on fc R,evnoldsHl99i|) . Except for special 
cases, the macroscopic equations that describe the momentum distributions and jump conditions 
cannot be solved analytically. 

Suppose particle acceleration is viewed microscopically and the individual particles are allowed 
to scatter kinematically where momentum and energy can be specifically tracked as the particle 
interacts with the shock environment. This proces s lends i tself nicely to computer techniques, 
specifically to a Monte Carlo technique developed bv lEllisonI lll98ll) over the last twenty years and 
is now a sophisticated tool for analyzing both nonrelativistic and relativistic shocks. 

The technique is briefly described here and will be discussed in more detail in a later chapter. 
The shock region is divided into a number of grid zones, many more than the few shown in Figure 
11.51 The model allows a number of particles to be injected into the shock environment far upstream, 
which simulates unshocked particles in the interstellar medium being subjected to a shock front. As 
the particles pass through the grid zones, the momentum and energy of the particles are tabulated 
as corresponding "fluxes" . In the case of so-called test particles, which do not modify the shock 
velocity profile, acceleration is turned off in the model, a compression ratio r — uq/u2 is chosen and 
the fluxes are tabulated at all the grid zones, both upstream and downstream. After the particles 
have completely passed through the shock and if the momentum and energy fluxes vary (i.e., are 
not conserved) across the shock, then the process is repeated with a different compression ratio until 
the fluxes appear to show continuity and remain at the same value upstream to downstream. 

The Monte Carlo model can also modify the shock velocity profile by allowing the backpressure 
of accelerated particles to affect the shock in a self-consistent manner. Again, after particles are 
injected, the Monte Carlo model measures the momentum and energy flux at each grid zone, com- 
pares the fluxes to those far upstream and estimates a new shock velocity at each grid zone that may 
better conserve momentum and energy flux. After a series of iterations, a final shock velocity profile 
will result that self-consistently conserves momentum and energy flux in the presence of accelerated 
particles, provided the correct compression ratio was chosen. If the fluxes do not balance across the 
shock, a new compression ratio is chosen and the process is repeated. Another feature of the model 
allows the possibility of accelerated particles to escape from the shock, carrying away momentum 
and energy flux, as discussed in iBerezhko fc Ellison (1999). These ideas will be discussed more fully 
in later chapters. 

1.5 Overview and Objectives 

There are two separate areas of research activity for this dissertation, and they are related in their 
common objective to understand and explain the observed radiation in gamma-ray bursts. 
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Figure 1.5: The shock envh'onmcnt showing the plasma regions divided up into grid zones for flux 
calculations by the Monte Carlo model. 
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The first research activity involves the ongoing deve lopment of a relativistic nonhn ear Monte 
Carlo model, based on an earlier nonrelativistic version ("Ellison. Barin g, fc Jonesl Il996|) . that will 
simulate the acceleration of charged particles by relativistic modified shocks with oblique magnetic 
fields when it is completed. It will be shown that in its current state of development, the model very 
satisfactorily simulates the collisionless scattering and acceleration of charged particles in relativistic 
modified shocks with the magnetic field parallel to the normal of the plane of the shock. The 
development of the model, its operation, and the parameter characteristics for unmodified shocks 
when shock speed, magnetic field angle and pressure anisotropy are varied, are discussed in Chapters 
2 through 5. Chapter 2 explains how and why the stress-energy tensor is used to handle momentum 
and energy fluxes. The fluid and electromagnetic stress-energy tensors are defined and a new equation 
of state is developed. The new relativistic magnetohydrodynamic flux relations for momentum and 
energy are derived in Chapter 3, along with the new jump conditions across the shock and a new 
method of estimating the adiabatic index at all shock speeds. Also included in Chapter 3 is a 
discussion of some of the compression ratio and magnetic field characteristics for unmodified shocks 
when the shock speed, magnetic field angle and pressure anisotropy are varied. Chapter 4 explains 
the new relativistic momentum transformations between reference frames and how modification 
of the shock velocity profile by energetic particles occurs. Chapter 5 describes the Monte Carlo 
technique for simulating the kinematic scattering and acceleration of particles across a shock front 
moving at relativistic speeds. An application to parallel, modified shocks is described and the 
resulting characteristics of the particle distributions, from nonrelativistic to highly relativistic speeds, 
are discussed. 

The second area of research activity utilizes the relativistic shock model to study lepton and 
baryon acceleration by relativistic modified shocks with the magnetic field parallel to the shock 
normal. Chapter 6 provides details of this study where the sensitivity of lepton injection and 
acceleration to subshock size and shock speed is explored to find the maximum possible energy 
efficiency of leptons. New results are presented, given the assumptions in this model and if particles 
are accelerated by diffusive shock acceleration, that suggest leptons can never carry enough energy 
to explain the observed gamma-ray burst spectra when the lepton to baryon number density is of 
the same order. Also, the lepton to baryon particle number density ratio is investigated to find 
the conditions under which leptons can carry the energy needed to produce the radiation observed 
in gamma-ray bursts. A new result is shown where the equipartition of energy between leptons 
and baryons is achieved in a relativistic shock when the lepton to baryon number density ratio is 
approximately 3 x 10^, independent of shock speed, at least in the range of shock Lorentz factors of 
2 to 12 used in this study. The results are discussed in the Conclusions in Chapter 7, along with a 
summary of all the accomplishments in the two research areas that include the development of the 
relativistic Monte Carlo model and the lepton-baryon acceleration study. 



Chapter 2 

The Stress-energy Tensor 



The stress-energy tensor describes the condition of a medium at any point in spacetime. One 
advantage to using a tensor formulation is that it is a ve ry con densed and convenient way to express 
physical laws. A second advantage, quoting TolmaJ l(l934l) . ". .the expression of a physical law 
by a tensor equation has exactly the same form in all coordinate systems..", due to the general 
transformation rules for tensors; i.e., 

T^r.-.^O (2.1) 
will be transformed into an equation of the same form 

= (2.2) 

when the spacetime coordinates are transformed from {x ^} to {x ^,x'^,x^,x 

The stress-energy tensor is composed of two primary parts; the fluid tensor and the electromag- 
netic tensor. Each tensor will be discussed in the sections below. Before discussing these tensors 
and the corresponding equation of state, some of their components require elaboration. 

2.1 The fluid tensor and the equation of state 

Aside from the first (T°°) component of the fluid tensor which is the total energy density in the 
proper frame, the other stress components of the fluid tensor consist of momenta divided by a unit 
area. If the normal of the unit area is in the direction of the momenta, the resulting component 
is called pressure. If the normal of the unit area is perpendicular to the momenta, the resulting 
component is a shear stress. 

Pressure is a macroscopic description of the particle momentum in a given region of space. If 
the pressure is a scalar quantity, it describes the average momentum or kinetic energy of an ideal 
fluid in the absence of electromagnetic fields. If the fluid is an ensemble of charged point particles 
(i.e., a plasma) in the presence of magnetic fields, the pressure may have different values in different 
directions. In this case, pressure is described by a tensor. In either case, an adiabatic compression or 
expansion of the region affects the pressure through an adiabatic index. A plasma is dominated by 
magnetic field effects because the particles move in such a way as to make electric fields negligible. 
Since a particle's energy does not ordinarily change by magnetic field defiections (assuming no 
radiation here), the adiabatic approximation is, in general, a good assumption. 
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2.1.1 The pressure tensor 



Consider a plasma with a magnetic field at some angle Qb with respect to the x axis in the xz plane. 
The pressure tensor will be constrained to the gyrotropic case; i.e., the stress can have one value 
parallel to the magnetic field vector , but can have a different value perpendicular to the magnetic 
field vector (with rotational symmetry about the magnetic field vector) as mentioned in equations 
(|B.7|) and IIB.8|l . Thus, we have the space tensor 



L J 7 



Pll 






















p 



(2.3) 



where subscript m refers to the magnetic axis coordinate system. Therefore, the pressure-stress 
tensor is diagonal using the magnetic axis, but rotation from the magnetic field coordinates to the 
xyz coordinates about the y ax is will produce a non-diagonal pressure-stress tensor as described by 
lEllison. Baring, fc Jonel l|l996>l . Their 3-dimensional pressure-stress tensor in the xyz plasma frame 
has components 



corresponding to 



P|| cos^Gs + Pl singes (Pl - P|| ) sin Gs cos 9 

P_i 

(Pj_ - PiOsinescoses 



P^icos^es + P|| singes 



where 
and 

The third equation completes the set 
P — P 




p ( cof? Qb - sin^ 8b i _ „ 



\ sin 9b cos 9 b 



2P^, cot (29s) 



(2.4) 

(2.5) 

(2.6) 
(2.7) 

(2.8) 



2.1.2 The fluid stress-energy tensor 

The fluid stress-energy tensor in the proper frame is defined with the following corresponding com- 
ponents: 



fluid 



1 











\ 
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^ XX 
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^ xz 
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^yy 
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p 

^ zx 





Pzz 


) 



(2.9) 



using Pxx, Pzz and P^z = Pzx from the previous section. 

^Note the small correction from the published reference. 
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The T'^" component, e, is th e total energy d ensity in the proper frame or plasma frame. The other 
components Pij are defined by 'Tolmai] l)l934(l as the "absolute stress" components in the proper 
frame. Pij is the force parallel to the i-axis exerted on a unit area normal to the j-axis. Hence, the 
diagonal components can be considered a pressure, but the off-axis components are shear stresses. 
The absolute stress components represent a different physical concept than the thermodynamic scalar 
pressure P. Scalar pressure P assumes a Maxwell-Boltzmann distribution and is Lorentz invariant. 
The non-thermal components of the pressure-stress tensor T"^ embodied in the fluid tensor above, in 
general, transform significantly and pick up momentum flux components in reference frames moving 
with respect to the proper frame. It may be noted that at nonrelativistic speeds, the pressure-stress 
tensor is invariant under Galilean transformations because force and area are invariant. 

A general adiabatic equation of state can be created from the conservation of energy density 
when oblique magnetic fields are present: 



3(r 



(2.10) 



where the trace of the pressure-stress tensor is Tr{T"'^ = P\\ + 2P±, T is the adiabatic index, and 
p is the rest mass density. Using equations (12. ()|) and (|2.7|) one can write 



TriV'^} _ 1 
3(r- 1) ^ 3(r- 1)' 



(P|| + 2Pj_ 



r - 1 



p.: 



^(2 cote -tan 6) 



(2.11) 



In terms of magnetic field, where cot Qb — it- and tan 6b = -j^ , this becomes 



1 



3(r- 1)^^11 ' ^^^^ ^ r-i 

Hence, the adiabatic equation of state, valid for both relativistic and nonrelativistic shocks, is 



2Pi 



1 



Prr 



PIE 

3 



Bz 
Br 



(2.12) 



1 



r - 1 



Px. 



p. 



xz 

3 



B, Br 



pc 



(2.13) 



2.1.3 Scalar pressure 

Scalar pressure is isotropic and Lorentz invariant. It is based on either the ideal gas law P — nkT 
where k is Boltzmann's constant, or it is calculated by averaging the squares of the individual particle 
momenta. When scalar pressure is presented as a three dimensional tensor, it is diagonal with equal 
components in all reference frames. The simplified equation of state, discussed in Appendix B, is 
stated again here: 

P . „2 



r- 1 



pc 



(2.14) 



The fluid tensor in the plasma frame becomes 



fluid 
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p J 



(2.15) 
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2.2 The Electromagnetic Stress-energy Tensor 

The electromagne tic field tensor composed by iLandau fc Lifshit3 lll962fl or the equivalent field- 
strength tensor of ljacksoni l)l975|) . each effectively using the Lorentz gauge, can be stated as 



and it's dual, 
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Ex 





~B, 


By 


Ey 


B, 





-Bx 




-By 


Bx 








—Bx 


-By 


-B, 


Bx 





E, 


-Ey 


By 


-E, 





Ex 




Ey 
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(2.16) 



(2.17) 



and from these field-strength tensors, the covariant form of the inhomogeneous Maxwell equations 
may be written as 

El"; = (2.18) 



and the homogeneous equations may be written as 



accompanied by the continuity equation 



•^,7 = 



(2.19) 
(2.20) 



where = {cpq, J) is the current four- vector, J is the current density three-space vector, and pq is 
the electric charge density. 

However, what is actually required is the electromagnetic energy momentum tensor or elec- 
tromagnetic stress-energy tensor. This tensor can be constru cted from the fi eld tensors above as, 
for example, IAppI fc Camenzindl ijlQSSf) did, but referring to iTolmanI lll984l) . the electromagnetic 
stress-energy tensor may be written directly as 



/ E^ 



-'em 



J_ 

47r 
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E X B 
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Qxx 


Qxy 


Qxz 


Qyx 


Qyy 


Qyz 


Q zx 


Q zy 


Qzz 



\ 



(2.21) 



where the Q's are the Maxwell stresses defined as: 

E^ +B^ 



Qi 



E 



Bf and Q,j = -{E.Ej + B.B^) 



(2.22) 



The overall electric field in a plasma is negligible and the plasma is dominated by magnetic fields. 
In addition, the xyz coordinate system is oriented such that the magnetic field lies in the xz plane; 
hence, the electromagnetic stress-energy tensor can be simplified to 



Chapter 2. The Stress-energy Tensor 15 
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-B^B, 
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(2.23) 





Bl - Bl 



2.3 Summary 

The concept of the stress-energy tensor was presented, and why and how it is used to handle 
momentum and energy fluxes, both fluid and electromagnetic, in the Monte Carlo relativistic shock 
model. The stress-energy tensor was shown to consist of two parts. The first part contains total 
fluid energy and a gyrotropic pressure tensor. The second part contains magnetic fields, with the 
assumption that electric fields, over the large scale, are negligible in the space plasma. Prom these 
tensors a new equation of state was established that is valid for both nonrelativistic and relativistic 
shocks. The equation of state includes a gyrotropic pressure tensor and oblique magnetic fields. 



Chapter 3 



Relativistic Magnet ohydro dynamic Jump 
Conditions 



Relativistic shock jump conditions have been presented in a variety of ways over the years. The 
standard technique for deriving the equations is to set the divergence of the stress-energy tensor 
equal to zero on a thin volume enclosing the shoc k plane and use Gauss's theorem to generate 
the jump conditions across the shock. For example, iTaubl l)l948|) developed the relativistic form of 
the Rankine-Hugoniot relations, using the stress-energy tensor with velocity expressed in ter ms of 
the Maxwell-Boltzmann distribution function for a simple gas. Ide Hoffmann fc pre- 
sented a relativistic MHD treatment of shocks in various orientations and a treatment of oblique 
shocks for the nonrelativistic case, eliminating the electric field by transforming to a frame where 
the flow velocity is parall el to the magnetic f i eld ve ctor (now called the de Hoffmann- Teller frame) . 
IPeacockl jrQSl V foUowinglLandau fc Lifshit3 l|l959l) . presented jump cond itions without elect r omag- 
neti c field s , and felandford fc McKeel l)l97fil) . also using the approach of iLandau fc Lifshitd l)l959l) 
andlTauhl ll 19481). develop e d a co ncise set of jump conditions for a simple gas using scalar pressure. 
IWebb. Zank fc McKenz"y l|l987j) provided a review of relativistic MHD shocks in ideal, perfectly 
conducting plasmas, and in particular the treatment bv iLichnerowic z' (1967. 1970*), which used this 
approach to develop the relat i vistic a nalog of Cabannes' s hock polar (|Cabannes,. .1970^1 whose origins 
also lie in lT.a,nda,u fc LifsTit^ llT95i. iKirk fc WebbI l(l988l) developed hydrodynamic equations using 
a pressure tensor, and IAppI fc Camenzindl l|l988l) developed relativistic shock equations for MHD 
jets using scalar pressure an d mag netic fields with components and (a parallel field with 
a twist). iBallard fc Heaven^ l)l99l|) derived MHD jump conditions using the stress-energy tensor 
with isotropic pressure and the Maxwell field tensor. By using a Lorentz transformation to the de 
Hoffman- Teller frame, they restricted shock speeds, wq, to uq/c < cosBbii, where Qbu is the angle 
between the shock normal and the upstream magnetic field; hence, this approach may only be used 
for mildly relativistic applications. All of these approaches assumed that particles encountered by 
the shock did not affect the shock velocity profile, i.e., shocked particles were treated as test particles. 

Here, previous work is extended by deriving a set of fully relativistic MHD jump conditions with 
gyrotropic pressure and oblique magnetic fields. The results are not restricted to the de Hoffmann- 
Teller frame and apply for arbitrary shock speeds and arbitrary shock obliquities. Solutions to the 
equations determine the downstream state of the gas in terms of the upstream state for the special 
case of isotropic pressure and, by parameterizing the pressure parallel and perpendicular to the 
magnetic field, for gyrotropic pressure. In this chapter analytic methods are used for determining 
the fluid and electromagnetic characteristics of the shock and do not inc lude particle acceleration . 
Later wor k will combine thes e results with Monte Carlo techniques (e.g., ison. Baring, fc JonesL 

that will allow the modeling of diffusive particle acceleration, including 
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the modification of the shock structure resuhing from the back-reaction of energetic particles on the 
upstream flow. 

3.1 Derivation of MHD jump conditions 
3.1.1 Steady-state, plane shock 

Utilizing a Cartesian coordinate system with the -l-x-axis pointing to the right, an infinite, steady- 
state, plane shock is shown in Figure IXTl traveling to the left at a speed uq with it's velocity vector 
parallel to the normal of the plane of the shock. In the rest frame of the shock, the upstream fluid 
appears to be flowing to the right with a speed of uq. 



Figure 3.1: Schematic diagram of a plane shock, in the rest frame of the shock, showing oblique 
magnetic fields and a Gaussian volume over which the divergence of the stress-energy tensor is 
integrated. In all of the examples the upstream flow is parallel to the shock normal. 

The upstream unshocked fluid consists of a tenuous, plasma of baryons and Icptons in thermal 
equilibrium with T^o = %iOj where Tpo (T^o) is the baryon (lepton) temperature. A uniform magnetic 
field, B, makes an angle Qbo with respect to the cc-axis as seen from the upstream plasma frame. The 
field is kept weak enough to insure high Alfven Mach numbers and thus to insure that the magnetic 
turbulence responsible for scattering the particles is frozen into the plasma. The xyz coordinate 
system is oriented such that there are only two components of magnetic field, B^o and B^o, in the 
upstream frame. The field will remain co-planar in t he downstre am frame and the downstream flow 
speed will be confined to the x-z plane as well fe.g.. I.Tones fc EUison. 1987) . 

In the shock frame, the upstream flow is in the -|-a;-direction and is described by the normalized 
velocity four-vector /3q = 70(1, /?q, 0, 0). The downstream (i.e., shocked) flow velocity four-vector is 
/32 — 72(1, 02, 0, /3|), where the subscript (2) refers, here and elsewhere, to upstream (downstream) 




shock front 



Chapter 3. Relativistic Magnetohydrodynamic Jump Conditions 



18 



quantities, P = u/c, and 7 = (1 — is the corresponding Lorentz factor associated with the 

magnitude of the flow velocity, u. 

The set of equations connecting the upstream and downstream regions of a shock consist of the 
continuity of particle number flux (for conserved particles), momentum and energy flux conserva- 
tion, plus electromagnetic boundary conditions at the shock interface, and the equation of state. 
The various parameters that define the state of the plasma, such as pressure and magnetic field, 
are determined in the plasma frame and must be Lorentz transformed to the shock frame where 
the jump conditions apply. In general, the six jump conditions plus the equation of state cannot be 
solved analytically because the adiabatic index (i.e., the ratio of specific heats in the nonrelativistic 
and ultrarelativistic limits) is a function of the downstre am plasma parameters, creating an inher- 
ently nonlinear problem fe.g.. l EUison fc Revnoldslll99l|) . Even with the assumption of gyrotropic 
pressure, there are more unknowns than there are equations; however, with additional assumptions, 
approximate analytic solutions may be obtained. 



3.1.2 Transformation Properties of the Stress-Energy Tensor 

Following iTolmanI l)l934|) . equations H2.9|l and (|2.23|) are combined to create the total stress-energy 



tensor, T^^tB-V sma of fluid and electromagnetic parts, i.e., 



^ total ~ ^ fluid 



(3.1) 



Then, by integrating over a thin volume containing the shock plane as shown in Figure 13.11 and, 
with Gauss' theorem, obtain the energy and momentum flux conditions across the plane of the 
shock by using T'^'^rii, = 0. Thus, it can be seen that T'^'^n„ yields the conservation of energy flux, 
T^^Ui, yields the a;-contribution to momentum flux conservation in the x-direction, T^^'rii, yields the 
z-contribution to momentum flux conservation in the x-direction, and rii, = (0, 1,0,0) is the unit 
four- vector along the x-axis in the reference frame of the shock. The Einstein summation convention 
is used here and throughout this thesis; Repeated Greek indices are summed over four-space(time) 
and repeated English indices are summed over three-space. 

The components of the fluid and electromagnetic tensors are defined in the local plasma frame 
and are Lorentz transformed to the shock frame where the flux conservation conditions apply, i.e.. 



1^ total is 



(3.2) 



where the subscript p (s) refers to the plasma (shock) frame. Since the flow speeds in the model 
may have two space components, in the x and 2-directions, the Lorentz transformation is 



A'' 



/ 7 




7 







7A 
7-1 







7-1 



^(izPx 1 + 



7-1 



PI 



(3.3) 



where 7 (1 - /32) - 2 , p^Pi is the square of the normalized flow speed as seen from the shock 
frame, and i is summed over the space components. 
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3.1.3 Flux Conservation Relations 

As discussed above, the energy and momentum conservation relations in the shock frame can be 
derived by applying T^'^rii, = to equation (|3.2|l . individually on the Lorcntz transformed fluid and 
electromagnetic tensors. 

The conservation of energy flux derives from 

{T^totai}, = {^dd}, + {T^m}, = ■ (3.4) 
The fluid contribution to energy (scalar) flux conservation is 

f!:'" = {Tfl^uid}, = j'PAe + P..) - 7(7 " l)^iP.. ' Pzz) + 



7[(27-l)/3^+/3.'] ^P.. , (3.5) 



while the electromagnetic contribution is 



Ff^ = {tO^}, ^^^^ [(7 - mPlBl + {iPl + PDP.BI- 

{[2-f~l)pl+l3l}PAB,] . (3.6) 

The conservation of momentum flux derives from 

{T^totai}, = {rZd}, + {^^m}, = . (3.7) 
The a;-component of the fluid tensor contributing to momentum flux is 

i^p"'' = {Tfiuid}, = I'Pli^ + P..) + P.. - (7 - 1?^{P.. P..) + 

2(7-l)(7/3^+/3f)^P.. , (3.8) 

and the z-component is 

F^f"" = {T^luid}, = l^PAie + P..) ^ {l - ml + lPl)^{P.. - P^z) + 

2 o2 



7 + 2(7-1)^^ 



Pxz . (3.9) 



The x-component of the electromagnetic tensor contributing to momentum flux is 
Ff^ = {T^^}, = ^PlB^ + ^ [{iPl + Plf - (7 - I? fulfil] {Bl - BI) - 



^(7 - l)(7/3^ + PD^B^B^ , (3.10) 



and the z-component is 

„2 



F^^ = {Ti^}, = i^P.P.B' + ^(7 - ifiPl - 0l)^{Bl Bl) 



i(7 - if^B^B, - ^B^B, . (3.11) 
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In all cases where the Alfven Mach number is large, the downstream flow velocity deviates only 
slightly from the shock normal direction so (3z ^ Px as shown in the bottom frame of Figure 13.21 
This allows a first-order approximation in Pz and the above equations become: 

^fluid ^ ^2^^(g ^ p^^) ^ ^(27 - l)l3zPxz , (3.12) 

F!"^ « ^PxB^z - 4^(27 - l)PABz , (3.13) 

F^^''' « I'Pli^ + P..) + P.. + 27(7 - 1)^^.. , (3.14) 

Ff^ - ^^PlB^ + iiBl - Bl) - ^A^^^B.Bz , (3.15) 



« l^PxPzie + Pxx) + iPxz - (7 - l)^{Pxx ~ p..) , (3.16) 



and 



Ff^ « ^P^P^B^ + ^(7 - iffjBl - Bl) ^B.Bz . (3.17) 

As shown in the bottom frame of Figure 13.31 the approximation Pz ^ Px becomes progressively 
better as the shock Lorentz factor increases. 

When the shock speed is ultrarelativistic, or when the magnetic field is parallel to the shock 
normal, Pz — Q and the above equations can be simplified further. 
For the ultrarelativistic case: 

Ff:^'' = l^Px{e + Pxx) , 



2 

pEM _ 7 a p2 
^ en ~ ^^rx-l^z ' 



and 



For parallel fields: 



and 



pEM = tI n2 

47r ^ 



^pz — ^^^x^z ■ 

F^^'^^j'PAe + Pxx) , 

pEM ^ r, 

en ^5 

Pp"'' = 7'/3'(e + P..) + P. 

pEM ^ r, 

pa: ; 

Pp^' = , 

;iEM 



3.18) 
3.19) 
3.20) 
3.21) 
3.22) 

3.23) 

3.24) 
3.25) 
3.26) 
3.27) 
3.28) 



Fi^/^ = (3.29) 
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3.1.4 Jump conditions 

The jump conditions consist of the energy and momentum flux conservation relations, the particle 
flux continuity, and the boundary conditions on the magnetic fleld. The conservation of particle 
number flux^ is 







where the brackets are used as an abbreviation for 



(3.30) 



(3.31) 



This jump condition, as well as the ones that follow, are written in the shock frame and, as always, 
the subscript (2) refers to upstream (downstream) quantities. The remaining jump conditions are: 



fluid 



EM 



■fluid 



F, 



EM 



and 



F 



fluid 



;iEM 



Adding the boundary conditions on the magnetic field, 



and 

completes the set of six jump conditions. 



1 2 



= 



(3.32) 
(3.33) 

(3.34) 

(3.35) 
(3.36) 



3.2 Solving the jump condition equations 



At this point there are eight unknown downstream quantities {(3x2 1 /3z2, Bx:2, Bz2, Pxx2, Pxz2, £2, 
and 77.2) and only six equations.^ Assuming isotropic pressure replaces Pxx and Pzz with P, sets 
Pxz — 0, and yields 

T^^L = ( e + P)u^n- + Pg^- + r^- , (3.37) 
where g^^ is the Minkowski metric fe.g.. iRvbicki fc LightmanLll979() 



/ -1 \ 
0100 
0010 

V 1 y 



(3.38) 



^Assuming no pair creation nor annihilation, which is reasonable for a steady-state model. 
^Pzz can be eliminated with equation 12.81 . 



Chapter 3. Relativistic Magnetohydrodynamic Jump Conditions 



22 



Unfortunately, this only removes one unknown and the equations are still under constrained. To 
proceed further, the equation of state is utilized to find an approximate expression for T. 

Oblique shock jump conditions cannot, in general, be solved analytically even for isotropic 
pressure because the downstream adiabatic index, depends on the total downstream energy 
density and the components of the pressure tensor (or scalar pressure), which are not known before 
the solution is obtained. The problem is inherently nonlinear except in the nonrelativistic and 
ultrarelativistic limits where r2 — 5/3 and 4/3, respectively. Furthermore, the gyrotropic pressure 
components are determined by the physics of the model and do not easily lend themselves to analytic 
interpretation, although Kirk & Webb (1988) provided equations based on a power law distribution 
in momentum for the pressure tensor components in the special case of a parallel relativistic shock 
with test particle diffusive shock acceleration. 

An excellent approximation can be obtained if Vth ^ uq, where Vth is the average thermal speed 
of the unshocked plasma. In this case, after scattering in the downstream frame, all particles have 

7rci = (1 - PLr'^^ , (3.39) 

where 

P.. - ^ = (3.40) 

C 1 - (JoP2 

is the relative (3 between the converging plasma frames. Using simple kinetic theory, it can be shown 
(see Appendix B) that the pressure 

P=^<p-v> (3.41) 

o 

where n is the particle number density, and p and v are the particle momentum and velocity, 
respectively. Then, with the approximation for particle velocity and the isotropic pressure version 
of the equation of state feauation l2.14|l . 



P (l/3)prcl^'rcl , , _ 7rcl/3i 



;2 
rcl 



e-pc^ (7rci-l)mc^ 3(7roi-l) 

or, 

r2 = . (3.43) 

37rcl 

This approximation allows a direct numerical solution for isotropic pressure, arbitrary obliquity, and 
arbitrary flow speed. It might be noted that equation H3.43|l provides an upper limit to the adiabatic 
index because any particles accelerated by the shock would tend to raise the average Lorentz factor 
of the particles and cause the adiabatic index to decrease slightly. 

In Figure 13.21 the results are shown as a function of /3o7o for a particular set of ambient pa- 
rameters as listed in the figure caption. The approximate equations (|3.12() - (|3.17|l have been used, 
although the exact equations can be solved if necessary. 

There are a number of important characteristics of the solution. First, the solution goes smoothly 
from fully nonrelativistic to ultrarelativistic shock speeds and obtains the canonical values for the 
compression ratio r = Pxo/Px2 = 4 for high Mach number, nonrelativistic shocks, and r — 3 for 
ultrarelativistic fiows. Three cases were shown with different upstream magnetic field obliquities: 
©BO = 0.1° (solid curves), 9bo = 10° (dashed curves), and 9bo = 60° (dotted curves). In all cases, 
the downstream magnetic field angle shifts toward Qb2 = 90° as the shock Lorentz factor increases, 
indicating the importance of treating oblique fields in highly relativistic shock acceleration. The 
value of r is weakly dependent on ©bo at nonrelativistic speeds, and essentially independent of ©bo 
(or other ambient parameters) at relativistic speeds. The compression ratio is within 10% of 3 for 



Chapter 3. Relativistic Magnetohydrodynamic Jump Conditions 



23 




Figure 3.2: Compression ratio, r, Qb2, 9u2, and r2 versus /3o7o- All examples have ng = 0.1 cm~'^, 
To = 10® K (with equal electron and proton temperatures), and Bq = 10 fiG. The solid curves have 
©BO = 0.1°, the dashed curves have Qbo — 10°, and the dotted curves have 6bo = 60°. In the 
bottom panel, the three cases are virtually identical as they are in the top panel for Po"fo > 0.02. 
T he solid dots in the bottom panel are values calculated using the Monte Carlo simulation described 
in lEllison fc Double! ll200"l . 
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7o > 2.1 and within 1% of 3 for 70 > 7. For the ambient conditions used in this example, the angle 
the downstream flow makes with the shock normal, Qu2, is small at all /3o7o (note the logarithmic 
scale for Qu2), consistent with the assumption that (3^ Px- 

In the bottom panel of Figure 13.21 r2 determined from equat ion 13.4311 is compared with a 
calculation using the Monte Carlo simulation described in detail in lEllison &: Double! l|2002fl (solid 
points). The Monte Carlo value for r2 is determined, without Fermi acceleration, directly from the 
shocked particle distribution using the r shown in the top panel. It is essentially identical to that 
given by equation (|3.43(l at all Lorentz factors confirming the utility of this approximation. As noted 
above, if Fermi acceleration does occur, r2 will approach 4/3 at lower /3o7o due to the contribution 
of energetic particles. 

For gyrotropic pressure, an additional assumption is needed to obtain an analytical solution and 
this assumption is taken to be 

P± = aP\\ , (3.44) 

where a is an arbitrary parameter and a = 1 gives isotropic pressure. Equation (|3.44ll allows the 
illustration of the effects of anisotropic pressure but is not suggested as a meaningful physical model. 
Using equations l|2.6;i . H2.7l and H3.44(l yields 

Pxx = P\\ (cos2 Gbh + a sin^ Obh) , (3.45) 

Pxz = P|| (a - 1) sin OBn cos QBn , (3.46) 

and the result is a closed set of equations for the jump conditions for shocks with gyrotropic pres- 
sure, arbitrary obliquity, and arbitrary flow speed. Results for various Obo's and a's are shown in 
Figure 13.31 (the other ambient parameters are the same as used for Figure I3.2|l . The solid curves 
(which are identical to the curves in the top panel of Figure have a = 1, the dashed curves have 
a = 0.7, and the dotted curves have a = 1.3. In all cases, the pressure in the unshocked, upstream 
plasma is taken to be isotropic and a 7^ 1 is only applied downstream. 

The effects of anisotropic pressure on the compression ratio come about mainly through changes 
to Pxx and this depends strongly on the downstream obliquity, Ob2- In the top panel of Figure 1^31 
©BO = 0.1° and 8b2 ^ 0.1° for nonrelativistic and mildly relativistic shock speeds. Therefore, at 
these speeds P^x ^ P\\ and since Py = P±/a, the fraction of downstream pressure in Pxx is inversely 
proportional to a. As shown in the Figure, r is less than the isotropic value, nso, for a = 0.7 and 
greater than for a — 1.3. As the shock becomes fully relativistic, however, 8b2 approaches 90° 
for any 8bo > (see the middle panel of Figure and Pxx ^ P±- When this is the case, the 
fraction of pressure in Pxx will be proportional to a and r < rjso for a > 1 and r > riso for a < 1. 

The transition where r crosses occurs at slower shock speeds as 8bo increases (middle panel 
of Figure until Qbo is large enough (bottom panel of Figure so no transition occurs. The 
change in Ob2 is relatively small for the examples shown in Figure but 0u2 changes significantly, 
as shown in the bottom panel for 8b2 = 10°. Despite the larger 8u2 for a = 0.7, the approximation 
Pz <C Px should still be valid. 

3.3 Correspondence with Nonrelativistic Jump Conditions 

The relativistic jump conditions [i.e., equations (|3.30l) - (|3.36() ] must correspond to the nonrelativistic 
jump conditions when the shock speed drops into the nonrelativistic range. Here, the pertinent 
relativistic equations are (|3.5(l - (|3.11(l and not their approximations because at nonrelativistic shock 
speeds, the z component of the downstream fiow is, in general, not negligible. Therefore, to reduce 
the relativistic jump conditions to their nonrelativistic counterparts, start by rewriting equations 
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Figure 3.3: The top three panels show the compression ratio, r, versus /3o7o for shocks with varying 
anisotropy, a, and obhquity, Qbo- The bottom panel shows the angle the downstream flow makes 
with the shock normal, 6u2, versus /3o7o. As in FigureESl "o = 0.1 cm"^. To = 10^ K (with equal 
electron and proton temperatures), and Bq = 10 /iG. The solid curves have a = 1, the dashed curves 
have a — 0.7, and the dotted curves have a = 1.3. 
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fl2 

(|3.5|) - (|3.11|) to second order in /3; i.e., 7 « 1 + Next, combine the resulting equations with the 

equation of state (|2.13|) and ehminatc the energy density e, then carefully watch the terms as /3 is 
reduced with prope r algebra to the nonrelat ivistic range. The resulting equations will match those 
published in .EUison. Baring, fc Jonej 1 199fil) . i.e. 



and^ 



1 



P +^ 



= 



PUxWz + Pxz 



BxBz 
Air 



Pxx^x ^" Pxz \ U 



Ux ( 2ga: 

3(r -i)\B, 



Bx 



13,1 2 , u^Bi 



UzBxBz 
An 



Qe 



(3.47) 
(3.48) 



(3.49) 



where Qesc is inserted to account for any energy flux lost due to escaping particles at a free escape 
boundary (FEB). The remaining jump conditions are straightforward. 



3.4 Summary 

The new relativistic magnetohydrodynamic flux relations for momentum and energy were derived, 
along with the new shock jump conditions. Setting the divergence of the stress-energy tensor equal 
to zero leads to the momentum and energy conservation laws across the shock. The equations were 
written in full and also to first order in small f3z and f3z — because /3z is usually small for relativistic 
shocks and f3z = for all parallel shocks. 

The continuity equation and electromagnetic boundary conditions, the new equation of state 
derived in the previous chapter, and a newly dervied approximation for the adiabatic index over the 
trans-relativistic range were combined with the conservation laws across the shock to create a set of 
general relativistic jump conditions relating the upstream and downstream regions across the plane 
of the shock. The equations were solved analytically for the case of isotropic pressure. For the case 
of gyrotropic pressure, an additional anisotropy parameter a was introduced to allow a solution. 

Some initial results were shown for unmodified shocks to demonstrate how the compression 
ratio, adiabatic index, downstream magnetic field angle and fiow angle vary with shock speed and 
upstream magnetic field. Also, the gyrotropic pressure tensor anisotropy was varied to show the 
effects on the compression ratio and field obliquity. 

Finally, a correspondence was established between the relativistic and nonrelativistic jump con- 
ditions. It was shown that all of the shock parameters in the relativistic shock model vary smoothly 
over the entire range of shock speeds from nonrelativistic to highly relativistic. This is a new feature, 
especially for modified shocks. 



Note a small typographical error in the energy flux equation in the referenced paper is corrected here. 



Chapter 4 



Particle Acceleration at Relativistic Shocks 



4.1 Comparisons between relativistic and nonrelativistic shocks 

Relativistic shocks, where the flow speed Lorentz factor 70 = [1 — (uo/c)^]~^/^ is significantly greater 
than 1, are likely to be much less common than nonrelativistic shocks, but may occur in extreme 
objects such as pulsar winds, hot spots in radio galaxies, and gamma-ray bursts (GRBs). Largely 
motivated by the appli cation to GRBs, relativistic shocks have received much attention by a num- 
ber of researchers (e.g.. teediiarz_^_Ostrowskil^ 19961 iKirk et all I2OO0I: ISchlickeiser fc Dermeilliooot 



lAchterberg et alll200l[ ^^ir'^atznCT^r^cK5^0^^ except for some preliminary work 

done over a decade ago ( Sch neider fc Kirk<1987tiEllisoilll991allJl and, aside from lEllison fc Double! 
(|25o2), current descriptions of relativistic shocks undergoing first-order Fermi acceleration are test 
particle analytical approxim ations that do not include the back-reaction of the accelerate d particles 
on the shock struc ture (e g., |PeacoclIll98lHHeavens fc DrurvLligSSHKirk fc Webbl . ll988|) : however, 
iKirk fc Schneider! ^ll987h^ and !Ballard fc Heavens! l!l99l!) used Monte Carlo techniques to calculate 
their test particle results. This may be a serious limitation of relativistic shock theory in applica- 
tions, such as GRBs, where high particle acceleration efficiencies are often assumed and test particles 
are, by their very definition, negligible. 

In coUisionless shocks, charged particles gain energy by scattering back and forth between the 
converging upstream and downstream plasmas. This basic physical process of diffusive or first- 
order Fermi shock acceleration, is the same in relativistic and nonrelativistic shocks. Differences in 
the mathematical description and outcome of the process occur, however, because energetic particle 
distributions are nearly isotropic in the shock reference frame in nonrelativistic shocks (where v ^ uo] 
V is t he individual particle speed), but are high ly anisotropic in relativistic shocks (since w ~ mo ~ c) 
(e.g.. !Peacocil!l98ltlKi7k fc SchneideR !l987b^ . 

The most important results from the theory of test-particle acceleration in ultrarelativistic shocks 
are: (i) regardless of the state of the unsho cked plasma, particles can pick up large amounts of energy 
AE 7q in their first shock crossing cycle l!Vietri!.!l995!). but will receive much srualler energy boosts 
{<Ef I Ei>^ 2) for subsequent crossing cycles (e.g.. !Gallant fc AchterberS !l999!: FAchterberg et al. , 
2001 ) ^; (ii) the s hock compression ratio, defined as r = uq/u2, tends to 3 as uq c (e.g.. !Peacock . 



198l!:!Kirll [ 1988), where U2 is the flow speed of the shocked plasma measured in the shock frame;^ 



^Ei (Ef) is the particle energy at the start (end) of an upstream to downstream to upstream (or a downstream to 
upstream to downstream) shock crossing cycle. 

■^Note that the density ratio across the relativistic shock P2/P0 7^ uo/u2, in contrast with nonrelativistic shocks, 
because the Lorentz factors associated with the relativistic flows modify the particle flux jump condition. Here and 
elsewhere the subscript (2) is used to indicate far upstream (downstream) values. 
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and (iii) a so-called 'universal' spectral index, a ^ 4.2 — 4.3 (in equation 15. 1|) exists in the limits of 
7o ^ 1 and 60 <^ 1, where 66 is the change in direction a particle's momentum vector makes at each 
pit ch angle scattering (e.p;.. iBednarz fc OstrowskL .1998 : Ac hterbcrg ct al.. 2001 ) . 

lEllison &: Doublg l)2002j) found that these results are modified in mildly relativistic shocks, even 
in the test-particle approximation, and in fully relativistic shocks (at least for 70 ^ 10) when the 
back-reaction of the accelerated particles is treated self-consistently, which causes the shock to 
smooth and the compression ratio to change from test-particle values. In mildly relativistic shocks, 
f{p) remains a power law in the test-particle approximation but both r and a depend on the shock 
Lorentz factor, 70. When efficient particle acceleration occurs in mildly relativistic shocks (i.e., 
7o <^ 3), large increases in r can result and a power law is no longer a good approximation to the 
spectral shape. In these cases, the compression ratio is determined by balancing the momentum 
and energy fluxes across the shock with the Monte Carlo simulation. For larger Lorentz factors, 
accelerated particles smooth the shock structure just as they do in slower shocks, but r approaches 
3 as 7o increases. In general, efficient particle acceleration results in spectra very different from the 
so-called 'universal' power law found in the test-particle approximation unless 70 ^ 10. 



4.2 Relativistic momentum transformations 

Elastic scattering and the subsequent changes in the direction^ of a particle's momentum takes place 
in the plasma frame, but all particle positions and momenta, and the corresponding jump conditions 
are handled in the shock frame; hence, relativistic frame transformations must be made from the 
shock frame to the plasma frame and from the plasma frame to the shock frame. 

Referring to Figure I^ITl the downstream plasma frame is moving with velocity U2 = xux2 + zuz2 
with respect to the shock frame, and a downstream particle has momentum p2 = xpx2 + yPy2 + zpz2 
in the downstream plasma frame, with the angle f between these two vectors, where x, y and z are 

unit vectors in the x, y and z directions. The particle momentum is defined as p — ^^mv where v is 

_ 1 

2 



the particle velocity, 7^ 



mc J 



1 



is the corresponding Lorentz factor, — pt+pt+P^ 



is the particle's rest mass, and the Lorentz factor of the relative frame velocity, "fu2 = ^1 — -f 

Before the relativistic transformations can be performed, the momentum components of the par- 
ticle parallel and perpendicular to the frame velocity vector U2 are required. The parallel momentum 

component is simply the scalar product u ■ p where the frame velocity unit vector is u = x-^ + z^, 

and it is understood that in these equations and the ones below where no number subscripts are 
used, all quantities are downstream plasma frame quantities. The momentum vector parallel to u is 

P|| = u(m • p) or 



P\\ 



+ z 



Px 



(4.1) 



The component of particle momentum perpendicular to the frame velocty vector is p± — p — p\\ or 













p±^ X 








'Px 













Pz 



ypy 



^The magnitude of momentum remains constant in an elastic collision. 
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Figure 4.1: Schematic diagram showing the shock frame and a downstream particle with momentum 
P2 in the downstream frame. The downstream frame is moving with velocity U2 with respect to the 
shock frame. 



1-1 If) 



(4.2) 



The relativistic momentum frame transformations from the downstream frame to the shock frame 
are: 

(4.3) 



Eu 



P\\,sk = lu 

and the invariant perpendicular component of momentum 

PA.,sk = P± ■ 

This leads to 



= lu [P\\ + ■yvmu] 



Psk = P\\,sk +P±.sk 



(4.4) 
(4.5) 



where 



P\\,sk = 7u 



, Px+ [ I Pz + Ivinu^, 



(4.6) 



The particle momentum, now seen from the shock frame, requires the corresponding x, y and z 
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components, where Px,sk — x ■ Psk, Py.sk — y 'Psk, and Pz.sk — z -psk- The resuhing components are: 



and 



Px.sk 



Pz,sk 



(7«-l)(^) +1 



Px + lulvmu^ + (7„ - 1) ( ) pz 



(7„-l)M^ +1 



Pz + lulv-muz + ilu-l) 



(4.7) 



(4.8) 



The y component of momentum will remain invariant under the Lorentz transformations. 

Finally, the momentum components are made dimensionless by dividing through by rripUo, where 
TOp is the proton rest mass and uo is the shock flow speed: 



Px.sk 



ilu - 1) 



Px 



Pz,sk 



(7u - 1) 



1 



Plsk 

Pz 



rripUo 



rUpUo 

^ Py 

nipUQ 

rUpUQ 



+ ilu - 1) 



«2 



Pz 



+ {lu - 1) 



~1F 



Px 



(4.9) 

(4.10) 
(4.11) 



The relativistic transformations from the shock frame to the downstream plasma frame are 
generated in the same way, using the transformation 



P\\ = 7« [P\\,k - 'lv,sk'musk] 



with the final normalized equations: 

2 



Px,pf 



(7n-l)(^) +1 



'yu'~fv,sk I / -1 \ / UxUz \ 

Px.sk —^UxirLpUo + (7„ - 1) — 2- pz^sk 



P^:Pf 



(7.-1) 



Py,pf 



Pz,sk 



, Py 
rUpUo 

lulv.skmUz 
TOpUo 



+ hu - 1) 



) Px.sk 



(4.12) 

(4.13) 

(4.14) 
(4.15) 



where it is understood that p'^ ^f^, p'y ^j, and p'^ [p'^ ^y, p'y and p'^ pj) are momentum components 
in the shock frame (plasma frame) normalized by rUpUQ, m is the mass of the particle under consider- 
ation, and jv,sk is the Lorentz factor associated with the particle velocity as seen in the shock frame. 
Equations for frame transformations between the upstream plasma frame and the shock frame are 
generated in the same manner. 



4.3 Method of calculating momentum and energy flux for 
parallel relativistic shocks 

lEllison fc Revnoldsl lll99lh describe the Monte Carlo procedure for numerically calculating the var- 
ious fluxes in the shock frame, showing the flux sums A(x,r) corresponding to the number flux. 
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B{x,r) corresponding to the momentum flux, and C{x,r) corresponding to the energy flux, where 
X refers to the x component of momentum measured in the shock frame, and r is the compression 
rati o U[)lu2- 

Ieii ison fc Revnoldsl lll99j) also describe the general process for flux conservation across the 
shock, assuming no explicit particle escape. A compression ratio r is chosen, the fluxes A,B, and 
C are calculated first for an unmodified shock, then the backpressure is used to modify the shock 
velocity profile through several iterations until the flux proflles converge to a stable solution. If the 
calculated fluxes are not constant across the shock, r is varied and the process is repeated until 
equations (|4.16lin7l and 14.181) are satisfied. 

Recalling the flux conservation relations for parallel shocks in Chapter 3, speciflcally, equations 
(|3.24|) - H3.29|l . where P^x is now called P and w = e + P, jump conditions are created across the 
shock. The pos sibility of escaping particle s which carry away particle, momentum, and energy fluxes 
as discussed in lBerezhko fc EUisorJ <)l999(l . is now included, to give the following equations: 

2 2 

l^wo^ +Po = llwx'^ +Px + Fp (4.16) 

Jo^qUq = ■^IwxUx + Fe (4-17) 

The t hird equation results from the requirement that the number of particles is conserved across the 
shock [Landau fc Lifshit j l)l959() ]: 



7onoUo = ixnxUx + (4-18) 

The three equations above represent the momentum density flux, the energy density flux, and 
the particle number density flux respectively, as viewed from the shock frame, while the number 
density n and the enthalpy w = e + P are both measured in their respective local plasma frames 
l)Ellison fc Revnoldslll99lf) . The terms Fp^ F^ and Fn represent escaping momentum, energy and 
particle fluxes, respectively, all referenced to the shock frame. Note that the Lorentz factors are 
referenced to the flow velocity (as a result of the three space components of velocity) . The subscript 
refers to the far upstream frame and the subscript x refers to any frame to the right of a grid 
plane, as shown in Figure [T31 Recall that the grid zones come from the original shock frame divisions 
along the x axis within which the flow velocity is constant. Also, if both sides of equation (|4.18() are 
multiplied by the particle rest mass the third equation could be written as 70^0*^0 = JxPxUx- 

In addition to the conservation relations, there is a relation that combines the adiabatic equation 
of state and the conservation of total energy, i.e. a parameterization of the ratio of specific heats 
appropriate to the region of interest, 

T=l + P/{e + pc^), (4.19) 

where P is the kinetic pressure, e is the total internal energy density including the rest mass energy, 
pc^ is the rest mass energy density, and T is the true ratio of specific heats under special conditions 
i Blandford fc McKeei Il973) . These ideas are discussed in more detail in Appendix B. Note that 



the compression ratio r = mo/m2 depends on F, making the solution to these equations nonlinear 
l)Ellison fc Revnoldd. Il99ll) . 

*One could say that -fvmp, where 7^ is associated with the particle speed, but viewed from the shock frame ■y^ is 
the same on both sides of the equation and drops out. 
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4.4 Method for modifying the relativistic shock velocity pro- 
file 

This section describes a method for modifying the shock velocity profile using the relativistic flux 
conservation relations, including the fluxes from escaping particles. 
From equations H4.f Gl and I4.f 7() . one can write 

7>.^+Px = i?(x,r), (4.20) 

and 

7Xwx = C(a;,r) (4.2f) 

where B(x,r) and C{x,r) are calculated as described bv lEllison fc Revnoldsl l)l99l[) . Enthalpy w 
can be eliminated by writing the last equation as 

= (4.22) 

and substituting this relation into equation 14.20(1 . Solving for pressure P we have 

P. = S(x,r) - 7^4^^^ = B{x,r) - C(x,r)^^ (4.23) 

where Ux{old) is the existing velocity at grid zone position x. 

So far, using Ux{old), we have an estimate of the pressure. Going back to equation (|4.16|l . it 
can be noted that the left hand side of this equation (far upstream from the shock) is a constant 
momentum flux because there are no escaping particles, so Kp is calculated: 

^^wo^ + Po ^ Kp. (4.24) 
c^ 

Then one can rewrite equation (|4.1()|) as 

^2 

Kp^-flwx^ + P,+Fp (4.25) 

The left hand side of equation ((4.17|l . far upstream where there are no escaping particles, yields a 
constant energy flux, so is calculated: 

j^WoUq = K,. (4.26) 

Then one can rewrite equation (|4.17() as 

Ke = ilw^Ua, + Fe (4.27) 

with 

= (4.28) 

Kp and K,, are numerically calculated in the same way as are B{x^ r) and C{x, r), but far upstream 
and far from all velocity variations where there are no escaping particles. Fp and Fg, the momentum 
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and energy fluxes of escaping particles, are numerically calculated separately from B{x,r) and 
C{x, r), but in the same way as the other values. 
Substituting H4.28|l into H4.25|l . we have 

u 

Kp^{K,-F,)^ + P, + Fp , (4.29) 
but is an estimate of the new flow velocity: 



Ux{new) . (4.30) 



This equation can now be combined with equation (|4.23(l to give the new flow velocity estimate 
in terms of the calculated flux values: 



Ux{new) 



Kp-Fp- B{x, r) + C{x, r) "\ ' 



(4.31) 



{K, - F, 

This last equation estimates the new velocity profile for the next iteration for a given compression 
ratio. After a number of iterations, the velocity and flux profiles should stop varying, i.e., Ux{new) 
should become closer and closer to the previously calculated Ux- If, after the series of iterations, 
the flux profiles are constant across the shock (i.e., from upstream to downstream), then the correct 
compression ratio was used. If the flux profiles are not constant, the compression ratio must be 
varied up or down. In general, if the downstream side shows a jump upwards in momentum flux 
after the iterations, it means the compression ratio was too high and needs to be lowered. 



4.5 Calculating F 

After the iterations are completed and, assuming conservation of the particle, momentum and energy 
fluxes has been achieved, the parameter T can be calculated. If the pressure an d energy density are 
known , equation (|4.19(l can be used to calculate F as shown in section 2.3 of lEllison fc Revnold^ 
For Lorentz factors of 10 or greater, the F derived in this manner should agree very closely 
with the F calculated from equation H3.43|l : 

r = 1 + /?2 = 1 + - (4.32) 



4.6 Summary 

Relativistic momentum transformation equations were developed that relate the upstream and down- 
stream plasma reference frames to the frame at rest with the shock. The equations allow oblique 
flows and arbitrary particle momentum directions. 

A method was developed for modifying the relativistic shock velocity profile by using the rela- 
tivistic flux conservation relations, with arbitrary magnetic field angle, to treat the momentum and 
energy fiux at every grid zone, plus the fluxes from escaping energetic particles. 



Chapter 5 



The Monte Carlo Technique and Computer 
Simulations 



The description of particle diffusion and energy gain is far more difficult when 70 3> 1 because the 
diffusion approximation, which requires nearly isotropic distribution functions, cannot be made. Be- 
cause of this, Monte Carlo simulations, where particle scattering and transport ar e treated explicitly , 
and which, in effect, solve the Boltzmann equation with collectiye scattering fe.g..lEllison fc Eichleil 
1^84; Kirk fc Schne ider. 1987b: EUiso n. Jones. & R,evnoldstll99(llEllison fc R,evnoldaLll99fHOstrowskil 
1991; .Bednarz fc OstrowskL ■1996.) . offer advantages over analytic methods. This is true in the 
test-particle approximation, where analytic results exist, but is even more important for nonlinear 
relativistic shocks. 

In nonrelativistic shock s, for v ^ un, a diffusion-convec t ion equation can be solved directly for 
infinite, plane shocks (e.s.. lAxford. Leer, fc SkadronI 1197^ iBlandford fc Ostriker. 1978'! . yielding 
the well-known result 

f{p)d^p(xp~'' d^p with cr = 3r/(r-l), (5.1) 

where r is the shock compression ratio, p is the momentum, and f{p) d^p is the number density 
of particles in d'^p. Equation H5.1|l is a steady-state, test-particle result with an undetermined 
normalization, but the spectral index, cr, in this limit is independent of the shock speed, mq, or any 
details of the scattering process as long as there is enough scattering to maintain isotropy in the 
local frame. To obtain an absolute injection efficiency, or to self-consistently describe the nonlinear 
back-reaction of accelerated particles on the shock structure (at least when the seed particles for 
acceleration are not fully relativistic to begin with), techniques which do not require v ^ uq must 
be used. Furthermore, for particles that do not obey w additional assumptions must be made 

for how these particles interact with the background magnetic wa ves and /or turbulence, i.e., the 
so-ca lled "injection problem" must be considered (see, for example, Ijones fc Ellis '^ ll 99ll: Im^W 
Il998() . The Monte Carlo techniques described here make the simple assumption that all particles, 
regardless of energy, interact in the same way, i.e., all particles scatter elastically and isotropically 
in the local plasma frame with a mean free path proportional to their gyroradius. These techniques 
and assumptions have been used to calculate nonlinear effects in nonrelativistic collisionless shocks 
for a number of years with good success comparing model results to spacecraft observations (e.g., 
lEllison fc Eichleil ll984l:lEllison. Mobhis. fc PaschmannLll99dlEllison. Jones. & Barind.ll999l) . 

Early work on relativistic shocks was mostly analytical in the tes t particle approximation (e.g., 
IBlandford fc McKeel Il97fil: IPeaco ck'. 'l98l'; 'Kirk & Schneider|, Il987at iHeavens fc Drurvi Il98fil) . al- 
though the analytical work of ISchneider fc Kirk (1987.1 explored modified shocks. Test-pa rticle 
Monte Carlo techniques for relativistic shocks were developed by iKirk fc Schneider! ( 1987b(> and 
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lEllison. Jones, fc Reynold^ l|l99Cl() for parallel, steady-state shocks, i.e., those where the shock nor- 
ma l is par allel to the upstream magnetic field, and extended to include oblique magnetic fields 
by lOstrowski (^1991,) . Some preliminary work on modified relativistic shocks using Monte Carlo 
techniques was done bv lEllisonllll991albl) . 

Monte Carlo techniques are used to model nonlinear particle acceleration in parallel collision- 
less shocks of various speeds, including mildly relativistic ones. When the acceleration is efiicient, 
the back-reaction of accelerated particles modifies the shock structure and causes the compression 
ratio, r, to increase above test-particle values. Modified shocks with small Lorentz factors can have 
compression ratios considerably greater than 3 and the momentum distribution of energetic parti- 
cles no longer follows a power law relation. These results may be important for the interpretation 
of gamma-ray bursts if mildly relativistic internal and/or afterglow shocks play an important role 
accelerating particles that produce the observed radiation. For 70 ^ 10, r approaches 3 and the so- 
called 'universal' test-particle result of N{E) cx E^'^-^ is obtained for sufficiently energetic particles. 
In all cases, the absolute normalization of the particle distribution follows directly from the model 
assumptions and is explicitly determined. 



5.1 Monte Carlo Model 

The t echni ques used here are essen t ially identical to those described in lEllison. Baring, fc JonesI 
l|l996|) and lEllison. Jones, fc Barind l)l999() . The differences are that the code has been made fully 
relativistic and only results for parallel shocks with pitch-angle diffusion are presented in this chapter. 
The code is steady-state, includes a uniform magnetic field, and moves particles in helical orbits, as 
shown in Figure I^TI The Alfven Mach number is assumed to be large, i.e., any effects from Alfven 
wave heating in the upstream precursor are neglected. This also means the second-order acceleration 
of particles scattering between oppositely propagating Alfven waves is neglected. Such an effect in 
relativisti c plasma s with strong magnetic fields is proposed fo r nonlinear particle acceleration in 

GRBs bv lPelle tier (1999') (see also IPelletier fc Marcowitbl Il998ll. 

The pitch angle diffusion is performed as described in ( Ellison. Jones, fc RevnoldS. Il99fl() and is 
shown in Figure That is, after a small increment of time, 8t, a particles' momentum vector, p, 
undergoes a small change in direction, 59. If the particle originally had a pitch angle, 9 (measured 
relative to the shock normal direction), it will have a new pitch angle 9' such that 



cos 9' = cos 9 cos 59 + \/l — cos^ 9 sin 56' cos , (5-2) 

where is the azimuth angle measured with respect to the original momentum direction. All angles 
are measured in the local plasma frame. 

If 59 is chosen randomly from a uniform distribution between and 59raax and is chosen from 
a uniform distribution between and 27r, the tip of the momentum vec tor, referring t o Figure 
will perform a random walk on a sphere of radius p. As shown bv .Ellison. Jones, fc Revnolda (1990,) . 
the angle 59max is determined by 

59^,^^ {Q5t/t,f' ^ {Un/Nf , (5.3) 

where N ~ Tgj 5t 3> 1 is the number of gyro-segments, dividing a gyro-period Tg = 2'Krg/v. The 
time tc is a "turn around" time defined as tc — X/ v, where A is the particle mean free path. The 
mean free path is taken to be proportional to the gyroradius Vg = pc/{QeB) (e is the electronic 
charge, Q is the ionic charge number, and B is the local uniform magnetic field), i.e., A — i^Vg, 
where 77 determines the strength of scattering. In all of the examples given here 77 = 1 is set, or, in 
other words, the strong scattering Bohm limit is assumed. 
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Figure 5.2: Conceptual diagram showing a particle with momentum p elastically scattering on 
magnetic turbulence through an angle SO. 
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For a downstream particle to return upstream, it's velocity vector must be directed within a cone 
with opening angle 62 such that \v2 cos02| > U2, where V2 and 82 are measured in the downstream 
frame and 62 — 0° is in the — x-direction, i.e., along the shock normal direction. For fully relativistic 
shocks with V2 — c and U2 — c/3, cos 02 <; 1/3 for a downstream particle to cross the shock into the 
upstream region. When the particle enters the upstream region it must satisfy essentially the same 
constraint, i.e., jwocos^ol > uq, where now vq and are measured in the upstream frame. "'^ Since 
both the particle and shock have high Lorentz factors, one can write 

(1-1/72)^/2 1-1/(27.') 

where jy = [1 ~ (vo/c)^]^^^^ is the particle Lorentz factor. Since cos^^o — 1 — for small 9q, we 
have 

(5.5) 



92 ^ L 

— 2 9 
7o Iv 



For ultrarelativistic particles with 7„ ^ 70, 6*0 — I/70 f^e.g.. ICallant fc Achterberel Il999l) . but 6*0 
can be much smaller for mildly relativistic particles. 

In order to re-cross into the downstream region, particles must scatter out of the upstream cone 
defined by 6*0 and Achterberg ct al. (2001.) show that most particles are only able to change the 
angle they make with the upstream directed shock normal by \S9\ ^ Oq before being sweep back 
downstream, making the distribution of shock crossing particles highly anisotropic. Therefore, if the 
shock Lorentz factor 70 ^ 1, a larger fraction of particles re-cross the shock into the downstream 
region with highly oblique angles (as measured in the shock frame) compared to lower speed shocks 
(see Figure lOI discussed below). Particles crossing at such obhque angl es receive smaller energy 
gains than would be the case for an isotropic pitch angle distribution and lAchterberg et alJ l|2001.) 
go on to show that <Ef/Ei>^ 2 for a shock crossing cycle (after the first one). 

With these concepts in mind, the requirement must be 50max < ^'o: or 

N > A^,„i„ = Unj^ . (5.6) 

The result (as shown in Figure lF3|l is that the power law spectral index, cr, asymptotically approaches 
a maximum value as N is increased. If N is less than the value required for convergence to the 
asympotic value (and the gyro-segments are too large) , the distribution will be flatter than produced 
with the convergent value of N because more particles are able to cross from upstream to downstream 
with 6'sk ^ 90° and receive unrealistically large energy boosts. This effect has long been known 
from the comparison of pitch-angle diffusion to larg e -angle scattering in relativistic shocks (e.g., 
iKirk fc Schneiderl Il987bt lEllison. Jones, fc Revnoldst Il99fi() . For all of the examples reported on 
here, A'^ is chosen large enough so it makes no difference if 60 is chosen uniformly between 0° and 
<J^max or if cos 69 is chosen uniformly between cosJ^max and 1. Figure [^31 shows how the results 
depend on N for shock speeds ranging from fully nonrelativistic to fully relativistic. In all cases, 
as N is increased the spectral index approaches a maximum and for 70 7 the well known result 
a = 4.2-4.3 is obtained. The fact that the computation time for the Monte Carlo simulation scales 
as N and N oc places limits on modeling ultrarelativistic shocks with this technique. 

In Figure 15.41 pitch-angle distributions (measured in the shock reference frame) of particles 
crossing the shock are compared. The curves are normalized such that the area under each curve 

^In the test-particle approximation, mo is just the shock speed. In nonlinear shocks, the flow speed just upstream 
from the subshock at 2: = will be less than the far upstream shock speed, mq, as measured in the shock reference 
frame. 
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Figure 5.3: Power law spectral index, a, versus number of gyro-segments, N, for various 
test-particle shocks as labeled. In all cases, as N is increased the spectral index converges asymp- 
totically and obtain the known results of cr ~ 4 for nonrelativistic shocks with r = 4, and a ~ 4.23 
for fully relativistic shocks with r = 3. The parameter A'min is defined in equation H5.6|l . 
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Figure 5.4: Distribution of the cosine of the pitch angle, i.e., Px/pt, for particles crossing a; = 0. 
The solid and dashed curves in the top panel are for a shock with Lorentz factor 70 = 2. The dotted 
curve in the top panel is for a nonrelativistic shock with speed uq = 5000 km s~^. The bottom panel 
shows curves for a fully rclativistic shock with 70 = 10. In all cases, the curves arc normalized so 
that the areas under them is 1 and the x-component of momentum, p^, is positive when directed 
downstream. The nonlinear (NL) and unmodified (UM) shock results are labeled. 
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equals one and the x-component of particle momentum in the shock frame, Px, is positive when 
directed downstream to the right (see Figure lB.lOl for the shock geometry). Here, pt is the magnitude 
of the total particle momentum also measured in the shock frame. In the top panel, unmodified 
(UM) test particle and nonlinear (NL) mildly relativistic shock (70 = 2) with a nonrelativistic one 
(uq = 5000 km s~^) are compared. Particles crossing the 70 = 2 shock are highly anisotropic with 
Px/Pt strongly peaked near 0.35. In the nonrelativistic shock, the particles are nearly isotropic 
except for a slight flux- weighting effect. There is little difference in the distributions between the 
UM and NL shocks. In the bottom panel, the pitch-angle distributions for UM and NL shocks with 
7o — 10 are shown. While the distributions arc somewhat more sharply peaked, they are quite 

similar to those for 70 = 2 and show relatively small variations between th e UM and NL shocks. 

The main difference between the present code and an earlier code used bv lElliso n. Jones. Sz Revnoldd 



(|199O0 to model test-particle relativistic shocks is that the previous code used a guiding center ap- 
proximation with an emphasis on large-angle scattering rather than the more explicit orbit calcu- 
lation of pitch-angle diffusion used here. Other than the far greater range in 70 and the nonlinear 
results presented here, the work of Ellison. Jones, & Reynolds (1990) is consistent with this work. 

The particle transport is performed as follows. Particles of some momentum, ppf, (measured 
in the local plasma frame) are injected far upstream from the shock and pitch-angle diffuse and 
convect until t hev cross a grid zone boundary , i.e., a plane with its normal parallel to the upstream 
flow speed (see lEllison. Baring. Sz Jonel . 1X9961 for a full discussion). For unmodified shocks, the flow 
speed changes from the upstream to downstream value in a discontinuous step at the shock boundary, 
but for nonlinear shocks the bulk flow speed changes in small steps, each separated by a grid zone 
boundary, from uq far upstream to 1*2 downstream. In the unmodified case, the shock thickness 
is essentially zero (i.e., shorter than the distance a particle diffuses in St) but in the nonlinear, 
modified case, the shock precursor extends over the entire region of varying bulk fiow speeds and 
a small scale "subshock" (at position a; = in this simulation) exists where most of the entropy 
production occurs. When a particle crosses a grid zone boundary, ppf is transformed to the new 
local frame moving with a new speed relative to the subshock and the particle continues to scatter 
and convect. Each particle is followed until it leaves the system in one of three ways. It can convect 
far downstream and not return to the subshock, it can obtain a momentum greater than some Pmax 
and be removed, or it can diffuse far enough upstream to cross an upstream free escape boundary 
(FEB ) and be removed. Bot h Pmax and the position of the FEB are free parameters in this model 
(see .Berezhko fc VolkL Il997t for a discussion of the self-consistent determination of the maximum 
particle energy in nonrelativistic supernova remnant shocks). For the nonlinear calculations, the 
shock structure and compression ratio are iterated until the number, momentum, and energy fluxes 
are conserv ed across the shock. Thi s procedure has been detailed many times for nonrelativistic 
shocks (see lEllison. Baring, fc Jone4 1996, and references th erein) and the modifications required 
for relativistic shocks were given in Tillison fc Revnoldal lll99ll) . 

To avoid excessive c omputation, a probability of return calculation as described in detail in 
lEUison. Jones, fc R,evnolds (,199(t() is used. That is, the standard expression obtained bv IPeacockl 
is used, 

Pn=('—^)\ (5.7) 

\V + U2 J 

to determine the probability, Pji, that a particle, having crossed a particular point in the uniform 
downstream flow, will return back across that point. Equation (|5.7|) is fully relativistic and inde- 
pendent of the diffusive properties of the particles as long as they are isotropic in the U2 frame. 
This isotropy is ensured by only applying equation H5.7|l once a particle has diffused several mean 
free paths downstream from the shock. For ultrarelativistic shocks with v ~ uq ~ c and r ~ 3, 
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equation H5.7|l gives Pr ~ 0.25. 

For clarity, it is noted that this is not the same probabihty, Vict, used by lAchterberg et alJ 



l|200H) to determine the test-particle power law index, i.e., 

where N{E) (x E"" and, for fully relativistic particles, s = a — 2. The quantity <Ef/Ei> is the 
a verage energy ratio for a particle undergoing a shock crossing cycle. Although not explicitly stated 
in lAchterberg et aTl l)200l[l . it is clear from the context that V^et is calculated just behind the shock 
where the particle distribution is highly anisotropic. In this case, particles that have just crossed 
from upstream to downstream will be more likely to recross back into the upstream region than 
indicated by equation H5.7|l because their pitch angles are more likely to be highly oblique relative 
to the shock normal than in the isotropic distributions further downstream (compare the solid or 
dashed curves to the d otted curve in the top panel of Figure ICTl discussed below). For 70 — 10, 
lAchterberg et alJ l)200l|) find Vrct = 0.435 ± 0.005 and <Ef/Ei>^ 1.97 ± 0.01 giving the standard 
result s = 2.230 ± 0.012. As shown in Figure j^J the immodified results are consistent with this 
spectral index for 70 ^ 7. 



5.2 Test-Particle Results 

In Figure [531 particle distributions for unmodified shocks are shown with speeds ranging from fully 
nonrelativistic {uq = 5000 km s~^) to mildly relativistic {f3o — uq/c = 0.5) to fully relativistic 
(70 = 10). The nonrelativistic distribution matches the standard test-particle Fermi result of cr = 4 
for r = 4 and the fully relativistic result is consistent with the well-known limit of cr ^ 4.2 — 4.3 
as 7o — > 00 for r = 3. In the trans-relativistic regime, both the compression ratio and the spectral 
index vary with uq. For the mq = 0.5c distribution shown in Figure 1531 the compression ratio was 
determined by balancing the mass, momentum, and energy fluxes across the shock in the test-particle 
li mit, i.e., by ignoring any eff ects from the accelerated particles. This technique is described in detail 
in lEllison fc Reynold j <|l99lh . It was found that r = 3.8 ± 0.1 and the spectral index is cr ~ 4.01, 
slightly flatter than 3r/(r — 1) ~ 4.07, the nonrelativistic result. 

Figure 15.61 shows the compression ratio as a function of /3o7o (solid curve) , still ignoring the 
effects of accelerated particles. The compression ratio is determined self-consistently by balancing 
the momentum and energy fluxes across the shock with no restriction on the shocked or unshocked 
adiabatic index. As expected, r decreases smoothly from 4 for fully nonrelativistic shocks to ~ 3 
for fully relativistic shocks. The power law index, cr, is also shown (dashed curve) and this varies 
slowly from cr = 4 to cr ~ 4.23 between the two extremes. For comparison, r (solid dots) is shown 
for shocks undergoing efficient particle acceleration. For these points, all shock parameters except 
/3o7o and Pmax are kept constant. For the nonlinear shocks, the maximum cutoff momentum is set to 
Pmax = 10'* /3o7o TT^pC in all cases. ^ The nonlinear results are discussed in detail below, but here only 
emphasize that r in nonlinear shocks will be larger than the test-particle value for 70 ^ 10 . The test- 
partic le res ults shown in F igure 1^?^ are in close agreement with recent analytic results of iKirk et all 
ll20n(t and lGallantI l)2002j) . It seems that no analytic results exist for nonlinear, relativistic shocks. 

In Figure l^n the average ratios of momenta are shown (measured in the local plasma frame) for 
particles executing upstream to downstream to upstream cycles across the shock, <pf /pi>u-d-m 

^The value of pmax can have a large influence on the shock characteristics at low /3o7o because r is large enough 
(and a small enough) that particles escaping at pmax are dynamically important. When /3o7o becomes large enough 
so that r drops below ~ 4, pmax is no longer an important parameter for the shock structure. 
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Figure 5.5: Particle spectra, p^fip), versus momentum for various immodificd shocks with speeds 
as indicated. The test-particle compression ratios, r, and spectral indices, a, are noted. The far 
upstream plasmas in the 5000 km s~^, and (3q = 0.5 shocks are thermal at 10^ K. The 70 = 10 shock 
has a far upstream plasma which is a (5-function at 1 MeV. All spectra are calculated at the shock, 
in the shock frame, and the relative normalization is arbitrary. 
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Figure 5.6: The solid line is the compression ratio, r, and the dashed line is the spectral index, 
a, for unmodified (i.e., test-particle) shocks. The solid dots show r for shocks undergoing efficient 
particle acceleration. In all cases, r is determined for each /3o7o by balancing the momentum and 
energy fluxes across the shock. The maximum cutoff momentum, Pmax, is unimportant for the 
test-particle shocks, but is set to 10^ /3o7o rnpC for the nonlinear shocks. 
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injection momentum. 
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and downstream to upstream to downstream cycles, <p//pi>d-u-d- These results are for 70 = 10 
(r = 3) and show a slight momentum dependence at low momenta but conve rge to <p f /pi>= 2.2±0.1 
at high momenta. This value is close to <Ef /Ei>= 1.97±0.01 reported bv lAchterberg et'a l. (2001) 
for 7o = 10. As mentioned above, in the first upstream to downstream to upstream cycle, particles 
achieve a large boost in momentum as indicated in the figure. 

The difference between this value of <Pf/pi> at large pi and that of lAchterberg et'al] l|200lh 
is greater than the uncertainties and probably stems from the different assumptions made in the 
simulations. As discussed above, <pf/pi> depends critically on the average angle a particle makes 
when crossing the shock. While the large majority of particles in these simulations gain energy 
by crossing from downstream to upstream and then immediately (within a few Sfs) re-crossing 
back into the downstream region at oblique angles, a few manage to diffuse farther upstream (see 
Figure below) . When these particles re-cross the shock into the downstream region, they can 
do so at flatter angles and receive larger energy gains. One might speculate that differences in how 
these few particles are treated in the simulations might produce the differences in <pf/pi>. 

As an illustration of how particles interact with nonrelativistic and relativistic unmodified shocks, 
trajectories for two individual particles are shown in Figures 15.81 and 15.91 The lower panel in each 
figure shows a trace of the particle trajectory and the upper panels show the particle momentum, 
always measured in the local plasma frame. For the nonrelativistic shock ^Figure l5.8|l . the speed 
of the particle is far greater than the shock speed and it diffuses easily on both sides of the shock. 
When it crosses x = 0, it does so nearly isotropically (except flux weighting makes crossings with 
flat trajectories slightly more likely) and essentially always gains momentum. The momentum gain 
in a single shock crossing is small, but a particle can stay in the system for many crossings. 

When the shock speed, wq, is close to c, the particle will be convected downstream much more 
rapidly than in nonrelativistic shocks and few particles will be able to cross the shock many times. 
However, downstream particles that do manage to cross the shock into the upstream region do so 
with much flatter trajectories, as discussed above, and can receive large momentum boosts in a 
single shock crossing due to the shock's Lorentz factor (note the logarithmic scale in the top panel 
of Figure 1^^. In a typical shock crossing cycle, downstream particles gain momentum when they 
cross into the upstream region (see positions labeled a and b in Figure RTIjl . lose momentum when 
they cross back downstream because they cross with oblique pitch angles, but end up with a net 
momentum boost. However, as shown by the position labeled c, it is possible for a particle to diffuse 
farther upstream before being convected back to the shock. In this case, it can cross the shock with 
a flat trajectory and gain momentum upon entering the downstream region. If the acceleration is 
efficient, the few particles that diffuse far upstream carry enough pressure to produce the shock 
smoothing discussed next. 

5.3 Non-Linear Results 

As noted bv lAchterberg et alJ (|200ltl . the large energy boost particles receive in their initial crossing 
of the shock provides a natural injection process for further acceleration and suggests that relativistic 
shocks may be efficient accelerators. However, just as with nonrelativistic shocks, efficient accelera- 
tion limits the use of test-particle approximations and requires that the nonlinear back-reaction of 
the accelerated particles be treated self-consistently (e.g., Jones & Ellison, 1991). These nonlinear 
effects will result in a smoothing of the shock and a change in the overall shock compression ratio, 
just as they do in nonrelativistic shocks. The differences between test-particle and nonlinear results, 
for the parameter ranges investigated here, are large enough to produce spectra noticeably different 
from the often quoted N{E) cx E^^-^ and to influence applications to GRB models where high shock 
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Figure 5.8: Particle trajectory (lower panel) and momentum (upper panel) in an unmodified shock 
of speed uq = 5000 km s~^. The momentum is calculated in the local plasma frame, either upstream 
or downstream from the shock. 
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Figure 5.9: Trajectory and momentuin for a particle in an unmodified sfiock with speed uq = 0.9c. 
Note that the horizontal axis is split at 1200 scatterings. Note also that the momentum scale is 
logarithmic because the momentum excursions are very large. 
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efficiencies are assumed. This is illustrated with two examples, one mildly relativistic (70 — 1.4) and 
one more fully relativistic (70 = 10). The far upstream conditions have relatively little influence on 
the results as long as 7p <C 70, where 7^ is the plasma frame Lorentz factor of the far upstream, 
injected particles. For concreteness, in both examples the far upstream plasma is taken to be a 
thermal distribution of leptons and baryons at a temperature of 10^ K, but only baryon acceleration 
is considered; the electrons only being included for charge neutrality. 

5.3.1 Mildly relativistic shock: 70 = 1.4 

Fiefure Pi. 101 shows unmodified and nonlinear shock structures for 70 — 1.4. The top panel shows 
7„(a;) u{x), where 7«(a:) = {1 — [u(a;)/c]^}~^/^, the middle panel is momentum flux, and the bottom 
panel is the energy flux, all scaled to far upstream values. All curves are plotted versus x, where 
a; = is the position of the sharp subshock. A logarithmic scale is used for x < — lOrg^o and a linear 
scale is used for x > —lOvg^o, where = mpUo/{eB).^ In each panel, the sohd curve is from an 
unmodified shock with r ~ 3.6 (see Figure l^?H| . while the dashed curve is the momentum and energy 
fiux conserving result. 

For the pitch angle diffusion model, where particles interact elastically and isotropically in the 
local frame according to equation H5.2|l and the discussion following it, particles are accelerated 
efficiently enough at the unmodified shock that the momentum and energy fiuxes are not conserved 
and rise well above the allowed far upstream values. In order to conserve these fluxes, the shock 
structure must be smoothed and the overall compression ratio increased above the test-particle 
value. The computational scheme calculates this compression ratio and flux conserving profile and 
the result is the dashed curve in the top panel with the corresponding momentum and energy fiuxes 
in the middle and bottom panels. 

The source of the non-conservation of momentum and energy is the efficient acceleration of 
particles by the sharp fiow speed discontinuity. While the actual injection and acceleration efficiency 
depends on the particular pitch angle diffusion model, once the scattering assumptions are made, the 
kinematics determine the injection and acceleration of the particles without additional parameters.^ 
Of course it would have been possible to make assumptions which resulted in an acceleration efficiency 
low enough that momentum and energy are approximately conserved without significantly smoothing 
the shock structure or changing the compression ratio from the test-particle value. For example, 
one could have only allowed shocked particles above some Lorentz factor 7i„j to recross into the 
upstream region or arbitrarily restricted the number of particles that recrossed into the upstream 
to a small fraction, /inj, of all downstream particles. By making /inj low enough or 7i,ij high enough 
one could make the efficiency as low as desired. However, there are at least three reasons for not 
making such an assumption. The first is that restricting the acceleration efficiency requires additional 
parameters (i.e., 7inj and/or /i„j) to those needed to describe pitch angle diffusion.^ The second is 
that models with inefficient acceleration will not help explain GRBs (or other objects) that require 
high efficiencies. If relativistic shocks are inefficient accelerators they are not very interesting. If 
they are efficient, they will have a qualitative resemblance to the results shown here even if, as is 
likely, the actual plasma processes are far more complex than the simple model we use. A third, 
perhaps less compelling, reason is that identical scattering assumptions as used here have been 

^In these parallel shock simulations where Alfven wave production is ignored, the magnetic field, B, has no other 
effect than setting arbitrary length and time scales. 

^Other input parameters, such as the Mach number and pmax or the position of the FEB, influence the accel- 
eration efficiency, but these are "environmental" parameters rather than parameters needed to describe the plasma 
interactions. 

^Alternatively, a far more complex model of the plasma interactions can be postulated than done here, inevitably 
requiring additional parameters (e.g.. lMalkovlll998ll . 
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Figure 5.10: Unmodified (UM) and nonlinear (NL) shock profiles, i.e., -fu{x)u{x), and momentum 
and energy fluxes versus position, x. All quantities are scaled to far upstream values and in all panels 
the solid curves are results from unmodified shocks and the dashed curves are nonlinear results. The 
NL momentum and energy fluxes arc 3 to 4% below the far upstream values because particles escape 
at a maximum momentum, Pmax = 10**/3o7o"i'pC ~ 10^^ eV/c. iV/Afmin ~ 10 have been used in both 
cases. 
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used for some time in nonrelativistic s hocks and shown to match both space craft observat i ons (e .g., 
lEUison. M obius. fc Paschmannl Il990l) and hybrid plasma simulations fe.g.. lEllison et all Il993(l of 
coUisionless shocks. 

The increase in compression ratio from r ~ 3.6 to ~ 4.5 shown in Figure 15. 101 comes about, in 
part, because the particles escaping at Pmax carry away momentum and energy fluxes which make 
the shocked plasma more compressible. This effect is countered to some degree by the fact that 
the self-consistent shock produces a downstream distribution with a smaller fraction of relativistic 
particles than the test-particle shock so that the downstream ratio of specific heats F > 4/3. This 
tends to produce a smaller compression ratio. The escaping fluxes show up as a lowering of the 
dashed curves below the far upstream values, as shown in the bottom two panels of Figure 15.101 
and amount to about 3% of the far upstream values for both momentum and energy. Including the 
escaping fluxes, the momentum and energy fluxes are conserved to within a few percent of the far 
upstream values. While the changes seen here are similar to those seen and d iscussed for many years 
in efficient, nonrelativistic shock acceleration fsee lBerezhko fc EllisonLri999l and references therein), 
it must be noted that there are no known analytic expressions relating escaping fluxes and F to r in 
this trans-relativistic regime. One obvious difference is that for nonrelativistic shocks, the escaping 
momentum flux is generally much less than the escaping energy flux (when both are meas ured as 
fractions of incoming flux) since pevl/{poUQ) ^ peWe/(PoWo) when Ve ^ uq (see lEUisonl 1 1 9 8 5l) . Here, 
Ve is the velocity of the escaping particle. In relativistic shocks, Ve ^ uq ^ c so the escaping fluxes 
are about equal as shown in Figure 15. 101 

In Figure \E.ll\ v'^ f(v) is plotted for the 70 ~ 1.4 shock. The solid curve is from an unmodified 
(UM) shock and the dashed curve is the nonlinear (NL) result. The shock smoothing and increase in 
r produce substantial differences in the spectra even though both shocks have exactly the same input 
conditions, (i) The overall normalization of the NL spectrum is less, reflecting the conservation of 
energy flux, (ii) The NL result has the distinctive concave curvature seen in nonrelativistic shocks 
stemming from the fact that higher momentum particles have a longer upstream diffusion length and 
get accelerated more efficiently than lower momentum particles in the smooth shock, (iii) The slope 
at the highest momentum in the NL spectrum reflects the overall compression ratio and is flatter 
than the TP spectrum because r is greater, (iv) The "thermal" peak is shifted to lower momentum 
in the NL result and contains a larger fraction of mildly relativistic particles than in the UM result, 
i.e., F ~ 1.41 for the NL shock compared to F ~ 1.36 for the UM shock. 

5.3.2 Fully relativistic, nonlinear shock: 70 = 10 

Figure 15. 121 shows results for 70 — 10. In all panels, the solid curves are for the unmodified shock, 
while the dashed curves are for the flux conserving smoothed shock, both with identical input 
parameters. As before, the acceleration is truncated with a Pmax = 10^/3o7oWpC ~ 9.3 x 10^^ eV/c. 
In Figure I5T^ u{x) and ju{x) are plotted separately. Even though the length-scale of the shock 
smoothing is only a few and the change in r = 3.03 ± 0.01 is small, they bring the momentum 
and energy fluxes into balance from values a factor of five too large in the unmodified shock. The less 
than 1% difference between r — 3.03±0.01 and the canonical value of 3 is significant; a self-consistent 
solution does not exist outside of this range. 

As discussed above, both the shock structure and the overall compression ratio may be modified 
when particle acceleration is efficient. For mildly relativistic shocks, such as the 70 = 1.4 example, 
a unique solution can be determined directly from the conserved momentum and energy fluxes. For 
larger 70 's however, the fluxes are less sensitive to changes in the structure and compression ratio 
and an additional constraint is needed to obtain a unique solution. This constraint is provided by 
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Figure 5.11: Particle distributions, p'^f{p), for the shocks shown in Figure IB. 101 The nonhnear 
spectrum (dashed curve) shows the distinctive concave shape seen in efficient nonrelativistic shock 
acceleration, and has a greater fraction of low momentum particles than the spectrum from the 
unmodified shock. As in Figure 15.51 the spectra are calculated at the shock in the shock frame 
and truncated with apmax- Unlike Figure 1531 the normalization here shows the actual acceleration 
efficiency. 
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Figure 5.12: Unmodified (UM: solid curves) and nonlinear (NL: dashed curves) shock profiles as 
in Figure for 70 = 10. The acceleration is truncated by Pmax — 10'*/3o7o'7ipC and N/Nn^in ~ 10 
in both cases. 
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Figure 5.13: The top panel is the flow speed at x normalized to the far upstream shock speed, uq 
versus x, for the 70 = 10 shocks shown in Fiefure [5.12l The bottom panel is the flow Lorentz factor, 
7„(a;), versus x. In both panels, the solid curves are the unmodified shock results with r = 3 and 
the dashed curves are the nonlinear results with r ~ 3.03 ± 0.01. 
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the relativistic jump conditions for momentum and energy, i.e. 



9 ? 

liowo^ + Po = ll2Wi-^ + P2 ; (5.9) 



7«o^oWo = 7«2^2U2 , (5.10) 

where w = e + P is the enthalpy density, e is the total energy density, P is the pressure, and it is 
assumed that escaping fluxes are negligible. The energy density and pressure are related through a 
combination of the adiabatic equation of state and the conservation of energy, i.e., 

P=(r-l)(e-pc2) , (5.11) 

where pc^ i s the rest mass energy density and F is, in special cases, the ratio of specific heats (e.g., 
lEllison fc^evnolds> .1991.1 . Dividing equation H5.9|) by equation (|5.1U|I yields 



7^o(eo + Pb)uo c2 ' 7^2(62 + P2)u2 



or, in terms of beta's, 



\eo + Po/ Po Ve2 + P2/ P2 



The second term on the left hand side is small compared to /3o for unshocked upstream particle 
temperatures less than lO^K and upstream particle densities less than 100 cm^^. Neglecting this 
term, equation H5.13|l becomes 

/^o = + f ^) ^ , (5.14) 

or, 



,62+^2. 

62/32 - /3o(e2 + i"2)/32 + ^^2 = , (5.15) 



which has the shock solution 



/?2 - /3o/r = ^ 

262 



/?o(e2+P2)-V^/3o' (62 + ^2) 



- 4eoPo 



(5.16) 



For 7o ^ 10, /3o — 1 and r ~ £2/ P2- Furthermore, if ultrarelativistic downstream particle speeds can 
be assumed so that 62 ^ P2C^ , equation (|5.16|) reduces to, 

^ (5.17) 



F2-I 



For F2 = 4/3,r~3, the standard ultrarelativistic result. 

To obtain a unique shock solution for 70 ^ 3, both the shock structure and the compression 
ratio, tmCi &re modified for the Monte Carlo simulation, P2 and 62 are calculated from the resultant 
downstream particle distributions, and tmc — r as determined from equation (|5.16() are checked. If 
''MC ^ the shock structure and tmc are varied until a consistent solution is found. For smaller 
7o's, equation H5.16|l doesn't apply because escaping fluxes become significant. In this case, however, 
the changes in the momentum and energy fluxes from changes in the shock structure and tmc a-^e 
large enough that a unique solution can be found easily, as in the 70 = 1.4 example. 
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Figure 5.14: Particle distributions, p'^f{p), for the shocks shown in Figures [5.121 and 15.131 with 
7o = 10. The spectra for the nonhnear (NL) and unmodified (UM) shocks are labeled and both are 
calculated at a; = in the shock frame. The light-weight solid lines show spectral indexes, a — 4.23. 
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Despite the fact that u{x) is modified on a fairly smaU length scale, the resultant particle 
distribution function is changed substantially, as indicated in Figure 15. 141 In this figure, the solid 
curve is from the unmodified shock and the dashed curve is from the nonlinear shock, both having 
exactly the same input conditions. The unmodified spectrum in Figure l5.14l is similar to that shown 
with a dashed curve in Figure 15.51 only now the far upstream plasma is taken to be a thermal 
gas at a temperature of 10^ K rather than a delta function distribution of particles with speeds, 
Dinj = (2i?inj/mp)^/^, with Einj = 1 MeV, as was assumed for the example in Figure 1^31 

The shock smoothing has caused the low energy portion of the distribution to steepen and the 
overall intensity to decrease, reflecting the fact that the NL spectrum conserves momentum and 
energy while the UM one doesn't. The peaks in the two distributions at low momenta are also very 
different, with the NL spectrum having a larger fraction of slower particles than the UM one. These 
peaks result from the first shock crossing where all particles receive a large energy gain. In the 
UM case, a far greater fraction of the accelerated downstream particles are able to receive further 
energization by recrossing back into the upstream region than in the NL shock. The different speed 
distributions result in different F's and it was found, by directly calculating F from the distributions, 
that FuM ^ 1-34 while Fnl ^ 1-36. 

5.3.3 Acceleration efficiency 

The absolute acceleration efficiency can be determined in the self-consistent, nonlinear examples 
directly from the particle distributions. In Figure 15.151 e(> p) is plotted, i.e., the fraction of kinetic 
energy density above a momentum p versus p for the 70 = 1.4 and 70 = 10 examples. The fraction 
of kinetic energy in the quasi-thermal part of the distribution can be determined from the relatively 
sharp fall off of the distributions at low momenta. This occurs near mpC for 70 — 1.4 and near 
10 mpC for 7o = 10. If the acceleration efficiency is somewhat arbitrarily defined for these two 
examples to be ei.4(> mpc) and eio(> lOrripc), respectively, the result is ei.4(> iripc) ~ 0.7 and 
eio(> lOrripc) ~ 0.6. Of course, the behavior of e(> p) depends on the particle spectrum from which 
it is derived and thus e(> p) is flatter for the 70 — 1.4 shock, with r ~ 4.5, than for the 70 — 10 
shock, where r ~ 3.03 ± 0.01. Furthermore, e(> p) depends strongly on Pmax in the 70 = 1.4 shock 
where particle escape plays an important role in determining r. The maximum momentum has little 
influence for 7 = 10 because of the steep spectrum. The fraction of kinetic energy density above 
lO'^mpC is about 0.3 for 70 — 1.4 and about 0.1 for 70 — 10. 

5.4 Summary 

Particles gain energy in coUisionless shocks by scattering nearly elastically off magnetic turbulence, 
back and forth between the converging plasmas upstream and downstream from the shock. While this 
basic shock acceleration physics is independent of the speed of the shock, the mathematical modeling 
of the process depends critically on whether or not the acceleration is efficient and whether or not 
particle speeds, v, are large compared to the shock speed, uq. The Monte Carlo techniques discussed 
here, which do not require v ^ uq, are well suited for the study of relativistic shocks, and for any 
shock where nonlinear effects are important and the energetic particles originate as thermal particles 
in the unshocked plasma. Except for computational limits, these techniques allow calculations of 
efficient particle acceleration in shocks of any Lorentz factor. 

As a check of the Monte Carlo code, it was demonstrated that the well known test-particle 
power laws can be obtained in fully nonrelativistic and ultrarelativistic shocks fFigures 15.31 and [5 . 5|l . 
In trans-relativistic shocks, however, no such canonical results exist because the shock compression 
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Figure 5.15: Acceleration efficiency, e(> p), defined as the fraction of total kinetic energy density 
above p versus p, is shown for two nonlinear cases. The sharp dropoff at low momenta indicates the 
extent of the 'thermal' peak. 
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ratio, r, depends on the upstream conditions in a nonli near fashion, even for test- particle shocks 
where all effects of accelerated particles are ignored (e.g., lEUison fc I^,evnoldsl. ll99l . it was shown 
how r, and the resulting power law spectral index, cr, vary through the trans-relativistic regime in 
Figure 1^21 where r has been determined by balancing the momentum and energy fluxes across the 
shock in an iterative process. These test-particle results with self-consistcntly determined compres- 
sion r atios are in close agreement with the recent analytic results of Kirk et aL (,2000) and Gallant 
l|2002h . 

Despite the fact that relativistic shock theory has concentrated almost exclusively on test-particle 
acceleration, it is likely that relativistic shocks are not test-particle but inject and accelerate particles 
efficiently. The reason is that regardless of the ambient far upstream conditions, particles that are 
overtaken by an ultrarelativistic shock will receive a large boost in energy ~ 7rei in their first shock 
crossing. Thus, virtually all of the particles in the downstream region of an unmodified shock are 
strongly relativistic with w ~ c as seen in the downstream frame. The ability to overtake the shock 
from downstream and be further accelerated depends only on the particle speed feauation l5.7|l and 
the presence of magnetic waves or turbulence with sufficient power in wavelengths on the order of 
the particle gyroradii to isotropize the downstream distributions. It is generally assumed that the 
necessary magnetic turbulence is self-generated and if enough turbulence is generated to scatter 
high momentum particles (with very low densities) that constitute a test-particle power law, there 
should be enough generated to isotropize lower momentum particles (which carry the bulk of the 
density). If acceleration can occur at all, it is believed that it is likely to occur efficiently making it 
necessary to calculate the shock structure and particle acceleration self-consistently. Furthermore, 
if relativistic shock theory is to be applied to gamma-ray bursts, where high conversion efficiencies 
are generally assumed, nonlinear effects must be calculated. 

When energetic particles are generated in sufficient numbers, the conservation of momentum 
and energy requires that their backpressure modify the shock structure. Two basic effects occur: 
a precursor is formed when the upstream plasma is slowed by the backpressure of the accelerated 
particles and the overall compression ratio changes from the test-particle value as a result of high 
energy particles escaping and/or a change in the shocked plasma's ratio of specific heats, F. As 
indicated by the 70 = 1.4 example (Section 5.3.1), mildly relativistic shocks act as nonrelativistic 
shocks, showing a dramatic weakening of the subshock combined with a large increase in r {Fig- 
ure l5.10|l . These changes result in a particle distribution which is both steeper than the test-particle 
power law at low momenta and flatter at high momenta ('Figure l5.11() . 

In faster shocks (i.e., 70 > 3), the initial test-particle spectrum is steep enough that particle 
escape is unimportant so only changes in F determine r feauation l5.16|l . In contrast to nonrelativistic 
shocks where the production of relativistic particles causes the compression ratio to increase, it is 
shown that r decreases smoothly to 3 as 70 increases and the fraction of fully relativistic shocked 
particles approaches one (Figure [^J. 

The most important result is that efficient, mildly relativistic shocks do not produce particle 
spectra close to the so-called 'universal' power law having a ~ 4.3. This may be important for 
the interpretation of gamma-ray bursts since the internal shocks assumed responsible for converting 
the bulk kinetic energy of the fireball into internal particle energy may be mildly relativistic and 
the external shocks, believed responsible for producing gamma-ray burst afterglows, will always 
go through a mildly relativistic phase (see IPiranL Il999l for a comprehensive review of gamma-ray 
bursts). 



Chapter 6 



Characteristics of lepton acceleration in relativistic 
parallel shocks 



The study of lepton (primarily electron and positron) acceleration in the vicinity of relativistic 
parallel shocks is part of a larger effort to investigate particle acceleration by relativistic oblique 
shocks using the equations and techniques discussed in earlier chapters. Here, the magnetic field is 
constrained to lie along the shock normal (a parallel shock), but shock modification and smoothing 
of the shock velocity profile by the backpressure of energetic particles is allowed. 

As mentioned in Chapter 1, it is believed that gamma-ray bursts are dominated by leptons 
through pair production in the early stages of the burst, with some unknown mix of baryons which 
is likely to be at a comparatively low density internal to the forward shock, but possibly at a higher 
density external to the shock. The gamma-ray burst spectra observed by earth satellites are generally 
believed to be created by leptons accelerated to high energies by the internal shocks^ for the initial 
intensity variations, and b y the forward shock to create the afterglow. 

Numerous papers re.g.. lBelAir978bHBlandford fc Eichle^llQsilLe^insonlllQQfiHSchmitz fc Chann 

I20n2^ have discussed the electron "injection problem" in which there may be problems in accelerating 
thermal electrons (pre-acceleration) to an energy where they can be further accelerated by the main 
shock. Electrons and baryons, assuming equal temperatures, have much different gyroradii resulting 
from the large mass difference between the two particles. The result is that baryons can be effec- 
tively scattered by the turbulence generated by the shock and self-generated Alfven waves because 
the wavelengths are of the same order as the baryon gyroradii. It is not clear how electrons interact 
with the shock, but one assumption that might be made is that electrons, with their much smaller 
gyroradii, undergo little scattering until they somehow increase the ir momeritum e nough to have 
gyroradii on the same order as the Alfven waves created by baryons ( Levinsonl . Il992l) . Recall, from 



sect ion 1.2 and also Chapter 4 , that elastic scattering is essential to first order Fermi acceleration. 

lEllison fc Revnoldsl l(l99ll) assume that Alfven waves with sufficiently short wavelengths exist 
and explore electron acceleration with various injection energ i es, but a number of author s (e.g., 
iGleevei Il984l iMcClements. et.al Il997t ISchlickeiser fc Dermeii 1200(11 ISchmitz fc Chanmanl [2002) 
have suggested other energy boosting mechanisms as a means of pre-accelerating electrons. For 
example, they have suggested various types of resonances with the magnetic field that would increase 
the kinetic energy of the thermal electr ons to the point where th ey c an continue to accelerate by 
the first order Fermi or diffusive process. iDieckmann. et.aP l)200Clf) and lSchmitz fc ChapmanI l)2002j) 
sug gest the Bunem an instability as an energy transfer mechanism. 

lElhsonl l ll98lh and others have proposed that a discontinuous collisionless unmodified subshock 

^ Refer to Figure lOl 
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exists along with the larger smoothed shock, and the subshock continuously accelerates the lower mo- 
mentum electrons to a point where they can interact more effectively with the smoothed shock. For 
relativistic shocks, all particles will have roughly an average Lorentz factor of 7ici and a momentum 
P — 7rei"^i'rci after the first shock crossing, as discussed in section 3.2. Hence, downstream electrons, 
after the first crossing, are more likely to have sufficiently long diffusion lengths and may then begin 
to participate more effectively in the Fermi acceleration process, in particular, the subshock. 

The diffusion length can be explained in the following way. Consider first a mean free path, 
A = ryrg, for coUisionless shocks, defining A as the average distance required to deflect a charged 
particle through 90 degrees by a magnetic field. Bohm diffusion (for the highest magnetic turbulence 
and the strongest scattering) is assumed when rj — 1. rg — pc/{qB) is the gyroradius, where p is the 
momentum of the particle, q is the unit electric charge, and B is the magnetic field strength. The 
diffusion coefficient, k = \v/2> (a choice; for example, see Jones & Ellison, 1991,) . is analogous to 
the standard diffusion coefficient used in collisional plasmas. Given these parameters, the diffusion 
length, L — k/wq cx pv, provides the overall length scale and it is roughly proportional to particle 
energy. 

It should be emphasized that the actual plasma processes in the vicinity of the shock are un- 
known and the resulting particle spectrum is dependent on the Monte Carlo model of the shock. 
A current limitation to the Monte Carlo shock model and the results on lepton characteristics pre- 
sented here, is the absence of a cross-shock potenti al in combination with a transverse magnetic field 
llZel'dovich fc Raized. Il967t I Jones fc EllisonLll987|) caused by the different mobilities of electrons and 
baryons as the particles cross the shock. In the absence of external electric fields, the potential dif- 
ference across the shock is determined by the kinetic energy of the electrons, therefore the extent of 
the departure from electrical neutrality is on the order of a Debye length. The cross-shock potential 
might have some infiuence on the dynamics of low-energy particles and the distribution of energy 
between leptons and baryons. Despite this, the Monte Carlo model, which explicitly simulates the 
particle kinematics, has been very successful over the past two decades in explaining observed spec- 
tra in nonrelativistic shocks. The more general relati vistic version of th is model, discussed in detail 
in Chapter 5, is in agreement with analytic results ijEUison fc Doublfil . 120021) and should, at least, 
provide useful information on the broad features of lepton spectral characteristics. 

Electron kinematics in coUisionless shocks are still poorly understood, and the purpose here is not 
to address details of the electron kinematics directly, but to understand the general characteristics 
of electron spectra, given shock speeds, subshock sizes and lepton densities as input parameters in 
the Monte Carlo program that models the particle dynamics. Although the results presented here 
are a small part of the larger study, it is hoped that they will contribute toward the ultimate goal 
of creating computer generated gamma-ray spectra that can be compared to actual observations. In 
this way, bounds may be placed on the gamma-ray burst model and increase our understanding of 
the physics of gamrna-ray bursts. 

IJones fc E1fisoTil (|l99lh note that there is considerable evidence that shocks can directly acceler- 
ate ambient ther mal p articles and, for this reason, the assumptions made here are the same as those 
made bv lEllisonI ljl98lh . It is assumed that in all shock fronts, some discontinuous subshock can exist 
which h as a small enoug h width to allow thermal electrons to accelerate by the first-order Fermi 
process. iTidman fc Kra ll ( 1971) propose that very fast shocks have a very narrow width, smaller 
than the diffusion length of thermal electrons, which suggests that, even after the shock is modified 
and smoothed by the backpressure of accelerated particles, thermal electrons see at least some small 
subshock embedded in the main shock (shown in Figure WA^ . and that some non- negligible fraction 
of electrons should be capable of scattering downstream and return to the (sub)shock for additional 
acceleration boosts until electrons have a large enough diffusion length to interact with the entire 
smoothed shock. 
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In this chapter, lepton and baryon momentum flux distributions will be compared as a function 
of the subshock fraction and as a function of shock speed. It will be shown how sensitive the lepton 
momentum distribution is to subshock size and to shock speeds over a range of mildly relativistic 
shocks. Furthermore, leptons and baryons will be compared as a function of lepton to baryon particle 
number density in terms of injection efhciency and energy efficiency. The equipartition density levels, 
where baryons and leptons share energy equally, will be determined for leptons of various masses, 
and it will be shown how a comparatively large particle density of leptons affect shock modification, 
as opposed to shock modification by baryons. 

6.1 Shock smoothing and the subshock 

Nonlinear interaction of accelerated energetic particle s on the shoc k velocity profile is well known 
(e.g.. Bell, 1987; Drurv. .198,1 .Jones fc EUison . 1991: Berezhko fc EU ison. 1999). In Figure lO an 
idealized modified coUisionless shock is shown in the shock reference frame. Far upstream, the flow 
velocity is the original shock velocity, uq. Energetic particles have diffusion lengths proportional 
to their momentum so the few most energetic particles, carrying significant momentum, create a 
backpressure that modifies the shock on a larger scale. The less energetic particles are more numerous 
and modify the shock on a smaller scale. In this way the precursor is formed as shown, down to 
ui. Below Ml, unaccelerated (thermal) particles have much smaller energies and they interact with 
the shock on a much smaller scale. The portion of the shock front between ui and U2 is called 
the subshock; i.e., that portion of the shock that has maintained approximately its original narrow 
width. 

Slower shocks are modified on a larger scale, as shown in Figure lfT^ because the energetic, highly 
relativistic particles have speeds much greater than the shock speed and many particles penetrate 
far into the upstream region, causing significant backpressure there. As the shock speed becomes 
more and more relativistic, downstream energetic particles scattering back towards the upstream are 
unlikely to cross the shock into the upstream region unless they lie within a critical angle 9q ^ I/70, 
as discussed in detail in section 5.1. Particles with trajectories at larger angles are scattered back 
downstream without crossing the shock (Refer to Figures 15.81 and I5.9|l . Hence, fewer energetic 
particles cross the shock from downstream to upstream, the backpressure from energetic particles 
is lower and the highly relativistic shocks are modified on a smaller scale, as shown in Figure 16.21 
as the dotted line. As a result, the modified velocity profiles of highly relativistic shocks appear to 
be much sharper than the velocity profiles of mildly relativistic shocks when viewed with the same 
distance scale. 

The smoothing of the shock in the precursor region is such that the momentum and energy 
flux are conserved across the shock as shown in Figure 16.31 This plot was created by allowing a 
large number of particles to run through the Monte Carlo model with a shock speed Lorentz factor 
of 5 and compression ratio of 3.1. The model ran through a number of iterations, and after each 
iteration the model compared the flux on both sides of the shock (actually on both sides of every grid 
zone, both upstream and downstream) where the flux differences were noted and used to slightly 
adjust the shock velocity profile. If the correct compression ratio is used, the momentum and energy 
flux differences will be negligible after the last iteration and the fluxes will appear flat over the 
entire range of x values, shown as the solid lines in the second and third frames of Figure The 
dashed lines, representing the first iteration and the unmodified shock, show how bad the mismatch 
is in momentum and energy flux across the shock. Momentum and energy are not conserved in an 
unmodified shock when real particles are accelerated. Details of the shock modification method can 
be found in section 4.3. 
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Figure 6.1: Idealized modified shock velocity profile with the main shock having compression ratio 
r — uq/u2, showing the precm'sor formed from the backreaction of accelerated particles, and the 
unmodified subshock with compression ratio — ui/u2. The velocities uo, ui, and U2 are all 
normalized to uo, and flow from left to right in the diagram. 
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Figure 6.2: Two modified shock velocity profiles with Lorentz factors of 1.25 and 8.0 as seen from 
the shock frame. The profile sharpens as the shock becomes more relativistic. The subshock is at wi . 

The downstream flow velocity, U2, is determined by the compression ratio which varies, depending 
on the speed of the shock. The shock is always modified in such a way as to conserve momentum 
and energy flux across the shock. 
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Figure 6.3: The shock velocity profile for a modified shock (solid lines) with no subshock, and 

Lorentz factor 70 = 5, plus the corresponding momentum and energy flux profiles displayed on loga- 
rithmic scales. These may be compared with momentum and energy flux profiles for the unmodified 
shock (dashed) where the fluxes are not conserved across the shock. 
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6.2 Sensitivity of lepton momentum flux distributions to sub- 
shock size 

To determine the sensitivity of lepton injection and acceleration characteristics to the size of the 
subshock in the shock velocity profile for mildly relativistic shocks, subshocks of various sizes were 
explicitly constructed with an input parameter as fractions of the difference between the uq and 
U2 flow speeds. After the subshock sizes were chosen, the shocks were smoothed, either manually 
or with a combination of iterative smoothing, to insure that momentum and energy fluxes were 
conserved, as shown in Figure 1^31 

Manual smoothing means that a smoothing value was set as an input parameter. An algorithm 
in the Monte Carlo model uses this parameter to create exponential smoothing of the shock velocity 
profile into the upstream region. The model generates momentum and energy flux values at every grid 
zone based on the amount of smoothing introduced. After observing the resulting flux profiles and 
their deviation from a straight line across the shock (which means a deviation from the conservation of 
momentum and energy), the input smoothing value is adjusted and the model is rerun, one iteration 
at a time, until an acceptable flatness of the momentum and energy flux profiles is achieved, similar 
to the solid line in Figure The manual method is less efficient and sometimes not as accurate 
as the automatic iterative method that the Monte Carlo model uses to adjust the velocity profile 
at every grid zone. However, manual smoothing is necessary at high Lorentz factors because the 
flux differences between grid zones become too small and the automatic smoothing becomes less 
accurate. 

"Leptons" (i.e., particles chosen to have mass nie/nip = 0.1, 0.01 and 0.001) were used for 
this study. Electrons with their actual mass were not used due to computer run time limitations. 
Referring to Figure IKjI and the related text, the Monte Carlo model simulates pitch angle scattering 
by altering the momentum direction of a particle through a small random angle 69 (which depends 
on the particle's gyroradius) in a small increment of time, 6t. The result is that lighter particles 
with smaller momenta and smaller gyroradii require many more time steps to accomplish the same 
amount of scattering as heavier particles, roughly as the inverse ratio of the masses. Hence, computer 
runs using very small mass particles become excessively long, given the present CPU speeds that are 
available. Therefore, the strategy was to use three larger lepton masses to expeditiously generate 
data, and then extrapolate the results down to those for the mass of the electron. The results are 
straightforward and shown in Figure I^TI 



6.2.1 Lepton flux distributions 

Momentum or energy flux distributions provide insight into how the particles are interacting with 
the shock and undergoing Fermi acceleration. The momentum flux distribution is measured by 
keeping track of the particles that cross a chosen boundary with the momentum associated with each 
parti cle, similar to the flux measurement at every grid zone as described in (.Ellison fc R^evnoldj . 
Il99lh . But, rather than simply summing the total flux in the x direction, the momentum for the flux 
distribution is made to be omnidirectional with a cosine function and is sorted into bins depending 
on the magnitude of the momentum. The result is a plot of the momentum flux distribution, f{p), 
as a function of momentum, p. Figure 16.41 shows a momentum flux distribution of particles with 
three different masses accelerated by an unmodified shock. An unmodified shock implies that the 
subshock size is the same size as the main shock, and implies that the entire shock width is much 
smaller than the diffusion length of all particles in the system, therefore particles of all energies 
(determined by the input parameters) can be accelerated by the shock. Also, because all particles 
see the same entire shock width, the momentum fiuxes all have the same slope. 
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In the top frame of Figure l6.4l the momentum flux distributions f{p) cx p^"^ , i.e., equation H1.12|l . 
of accelerated baryons and leptons are shown simply as the logarithms of f{p) vs p. Their slope is 
called the spectral index, a. In the bottom frame, the distributions are flattened by multiplying them 
by p"^ , which allows more direct comparisons of subtle variations in the momentum distributions. 
The point here is to show the concept for unmodified shocks, and the flattening will become more 
useful for modified nonlinear shocks to be discussed later. 

As the particles scatter back and forth across the shock, they pick up increments of momentum, 
depending on their pitch angle at the time of crossing. Recalling the discussion at the beginning 
of this chapter, all particles pick up average energy of 7i.oimc^ and average momentum of 7rci"T'i^roi 
at the first shock crossing, hence the more massive particles initially have, on the average, more 
momentum than less massive particles in the ratio of their masses. For any given particle speed, the 
probability of return of the particle from downstream depends almost entirely on the compression 
ratio across the shock^ . Since all particles have the same probability of return for a given particle 
speed, the momentum distributions as a function of momentum will all have the same slope, as 
described in section 1.2. However, the distribution flux levels for the species were set at the earliest 
stages of acceleration, referring to equation Hl.ll|l . and they follow their corresponding flux levels 
for all higher momenta. For example, in the top frame of Figure the less massive leptons (i.e., 
me = 0.1 and 0.01) scatter and accelerate at lower momenta with initially large momentum flux 
[again, equation Hl.ll|l ] , but due to the constant spectral index, there are significantly fewer leptons 
and correspondingly lower lepton momentum flux at the initial momentum of accelerated baryons. 

Referring to Figure shocks modified by the backpressure of energetic baryons are smoothed 
and the size of the subshock is correspondingly reduced to ui. As a consequence, the less massive 
leptons, with initially less momentum and shorter diffusion lengths than baryons, cannot interact 
effectively with the entire shock because the shock is smoothed on the longer scale of the baryons. 
As a consequence, far fewer leptons reach energies where their diffusion lengths are comparable to 
those of baryons and the overall lepton flux is much lower. This effect can be observed in Figure lF31 
The leptons with mass me = 0.01 are affected much more drastically than the leptons with mass 
rrie —0.1 by the smoothed shock. 

To determine the sensitivity of the lepton momentum flux to the size of the subshock, subshocks 
of various sizes were explicitly constructed with an input parameter, as described above, in smoothed 
shocks with various Lorentz factors. The subshock size is deflned in the following way. Referring 
to Figure |0] again, the subshock fraction of the main shock is fs — (ui — U2)/{uo — M2). The 
compression ratio of the subshock is — Ui/u2 and the compression ratio of the main shock is 
r — uo/u2', hence, fs — (r-g — l)/(r — 1), and fs is the normalized or relative subshock size. Subshock 
size can also be stated simply by just comparing compression ratios, r^/r; however, this method 
gives a value of 1 when the subshock disappears at ui — U2 and does not show the fraction of the 
main shock ascribed to the subshock. 

The velocity profiles containing constructed subshocks were adjusted to conserve momentum 
and energy fluxes, using baryons only. The resulting velocity proflles were then used as a template 
to run the model with pairs of particle species, the first being baryons with mass equal to 1.0 and 
the second being lighter particles with mass of either 0.1, 0.01, or 0.001, each with its own number 
density of 1 particle per cubic centimeter. Little error was introduced by this method because the 
baryons are primarily responsible for the shock smoothing, and the smaller mass leptons were shown 
to have only a minor effect on the overall shock velocity profile smoothing and the related momentum 
and energy flux conservation proflles when the densities of baryons and leptons are approximately 
equal. 

^Note equation il.6i where U2 is determined by the compression ratio r. 
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Figure 6.4: Spectra for an unmodified relativistic (70 = 2) shock. The baryon and lepton spectra 
have the same slope, indicating that all particles have a sufficient diffusion length to interact with 
the shock. The bottom frame shows the same spectra, but flattened by multiplying the distribution 
function by normalized momentum raised to the power of a. 
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Figure 6.5: Spectra for the modified relativistic (70 = 5) shock shown in Figure The baryon 
and lepton spectra develop the same slope, indicating that the leptons have achieved a long enough 
diffusion length to interact with the entire modified shock, but at greatly reduced numbers compared 
to baryons. The bottom frame shows the same spectra, but fiattened by multiplying the distribution 
function by normalized momentum raised to the power of a. 
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After each run was completed, the flux spectra was generated and measured from plots similar 
to the one shown in the bottom frame of Figure 16.51 The vertical axes of the plots show the 
logarithms of the fluxes, therefore the differences of the measured values between baryons and the 
lighter leptons directly give the logarithm of the ratio of their flux distributions, logiQ[f{e)/f{p)]. 
A number of subshock sizes were explored, from = to = 0.8, for several shock speeds with 
Lorentz factors ranging from 70 = 2 to 70 = 12. The results, shown in Figure IfT^ indicate a great 
increase in sensitivity to subshock size as the particles become less massive, but show little sensitivity 
to shock speed for the larger subshock sizes over the range studied here. For large subshock size, 
the flux values approach those of an unmodified shock, such as those shown in Figure For small 
subshock size, lepton flux is extremely small and, as shown in Figure IFTfl with projections to the 
electron rest mass, rUf^/mp = 1836, electrons would require a significant subshock size, or would 
need to be initially pre-energized by some other means, in order to enter into the first-order Fermi 
shock acceleration process. However, highly relativistic shocks tend to become very sharp on the 
scale of baryon diffusion lengths, as shown in Figure 1^7^ and, for ultrarelativistic shocks, may be 
sharp on the scale of electron diffusion lengths and thereby have a large effective subshock capable of 
accelerating electrons. This may be the reason that there appears to be little dependence on shock 
speed when the constructed subshock is made large. 



6.2.2 Lepton injection efficiency 

The lepton momentum flux is directly related to the injection efficiency of l eptons taking part in the 
first-o rder Fermi acceleration process. Based on the concept expressed bv |EllisQn. Baring, fc Joiie^ 
l|l995|) . the relative injection efficiency, einj(p > Po)i can be defined as the number density of particles 
accelerated above a critical momentum, po, compared to the number density of particles no flowing 
toward the shock, or 

N _ n-mjiP > Po) (n IN 

emj (> Po) = — ■ (6.1) 

The critical momentum, shown as the dotted line in Figure I^TSI is approximately at p/nipC — 5 (i.e., 
logj^Q [p/mpc] « 0.7). This is the momentum where the thermal distribution of baryons ends and the 
power law distribution begins. The same momentum can be found in the lower frame of Figure 1^751 
at the beginning of the flat portion of the baryon flux distribution and where all three momentum 
flux distributions flatten to the same power law for this particular case of a modified relativistic 
shock with Lorentz factor of 70 — 5. The efficiencies of each type of particle was normalized to 1.0 
to allow the injection efficiency of each type to be measured against its own initial particle density. 
This also allows direct comparisons between injection efficiencies of the various particle types. 

The injection efficiencies for leptons and baryons may be determined as follows. By using the 
critical momentum, pq, as the defining point for baryons, the absolute injection efficiency for baryons 
is read off the plot shown in Figure I^TSl Then, using the same momentum, for leptons, the relative 
injection efficiencies for leptons are the differences (in logarithms) between the lepton efficiencies 
and the absolute baryon efficiency at pq. 

The critical momentum was found for each subshock size and for each Lorentz factor to provide 
the absolute injection efficiencies for baryons. Then, all of the relative injection efficiencies were 
found for leptons with masses of 0.1, 0.01, and 0.001 to create the injection efficiencies shown as a 
function of fractional subshock size, fs, in Figure The lepton efficiencies are normalized to the 
baryon injection efficiency, and the baryon efficiency is shown as the absolute efficiency on the same 
graph. The normalization factor is the ratio of the baryon to lepton densities, taken in pairs. The 
injection efficiencies tend to follow the same pattern as did the fluxes. This is expected because the 
fluxes above the critical momentum generally follow the same power law for all particles. It may 
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Figure 6.6: The ratio of lepton to baryon momentum flux as a fimction of subshock size is shown for 
particles with masses of 0.1, 0.01, and 0.001. The lighter leptons are far more sensitive to subshock 
size than the heavier lepton, spanning fifteen orders of magnitude for small subshocks. Sensitivity 
to shock speed is low except for very small subshocks where the flux is seen to increase with higher 
Lorentz factors. 
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Figure 6.7 : Projections of lepton flux from larger lepton masses down to the mass of the electron 
show what flux to expect for real electrons. Note that the subshock sizes are in terms of the 
compression ratios, rg/r. 



Chapter 6. Characteristics of lepton acceleration in relativistic parallel shocks 



72 



be noted in Figure 16.91 that the lepton injection efSciencies tend toward those of test particles in 
unmodified shocks (i.e., where fs — I) at the highest subshock sizes at all shock speeds used in this 
study. 

Also note, at the smallest subshock sizes, there is a greater dependency on shock speed. This 
characteristic is shown more clearly in Figure lB.lOl This figure shows the relative injection efhciency 
for leptons with masses O.f and 0.01 as a function of shock speed in shocks where there was no 
explicit subshock established. For the shock speeds with Lorentz factors of 2 - 20, there is a dramatic 
improvement in lepton injection efficiency through eight orders of magnitude in the me/nip = 0.01 
mass. Shock speed dependence occurs because, as the speed increases, the shock front becomes 
sharper on the length scale of the lepton and more of the lower energy leptons can enter into the 
acceleration process. Once a particle reaches a momentum where its diffusion length is large enough 
to interact with the entire smoothed shock width (i.e., through the entire precursor region in Figure 
16. 1|) . the power law portion of the flux distribution curve begins and the injection efficiency is 
established. Baryons reach this "critical momentum" point easily, compared to leptons. Lepton 
momentum is lower than baryon momentum at the same speeds due to their smaller masses. At 
shocks with low Lorentz factors, the diffusion lengths of all leptons increase substantially after the 
first crossing, but not enough to interact with the entire smoothed shock width. Leptons require 
additional shock crossings to reach the required energies to interact with the entire shock width 
and far fewer leptons can accomplish this, therefore the lepton injection efficiency is lower. The 
theoretical maximum relative injection efficiency would occur when the shock front is so narrow 
that the diffusion lengths of all particles, leptons included, is greater than the shock width after 
the first crossing. By the same argument presented earlier, the ratio of lepton to baryon momenta 
is lower by the ratio of their masses after the first crossing. Then, given equal number densities 
of leptons and baryons, the maximum relative injection efficiency of leptons depends only on the 
injection efficiency power law through the relation 

logio[einj(max)] = cri„jlogio^ . (6.2) 

I Up 

The average measured crinj ^ 1.3±0.1. Computer results for leptons with mass me/nip = 0.01 extend 
only to a Lorentz factor of 20, therefore it is not clear how their injection efficiency behaves at higher 
shock speeds, but based on the relationship above, the maximum relative injection efficiency for these 
particles may be approximately 10~^'^ at some Lorentz factor above 100 in Figure 1^.101 It appears 
that the maximum relative injection efficiency (i.e., at the same momentum) of leptons will always 
be far below that of baryons when their number densities are equal. 

6.3 Characteristics of increased lepton densities 

Injection efficiencies, as defined in the previous section show the fraction of particles that are incor- 
porated into the Fermi acceleration process, but a more important question is how much energy is 
carried by the particles. Present gamma-ray burst models require leptons to carry a significant frac- 
tion of the total energy density of the burst, and there are a number of ways this cound occur. One 
possibility is through rapid energy transfer from baryons to leptons by some unknown means. An- 
other possibility is through a mechanism other than shocks to directly energize leptons; for example, 
magnetic reconnection. The possibility that seems most likely, and the one that will be investigated 
here, is that there are more leptons than baryons, and they are energized by diffusive shock acceler- 
ation. One of the requirements in the Fireball model of a gamma-ray b urst is that the baryon m ass 
fraction must be low to allow the fireball to expand relativistically fe-g.- lMeszaros fc Reesl Il997|) . In 
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Figure 6.8: Injection efficiency einj{> p), i-C, the fractional number density of particles above a 
given momentum, is plotted against normalized momentum. The critical momentum, po, is chosen 
at the point where baryons clearly obtain a power law dependence. Lepton power law dependence 
is also established around po- 
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Figure 6.9: The top curve shows the absolute injection cfRciency of baryons in modified shocks 
as a function of subshock size. The curves for rUe = 0.1, 0.01, and 0.001 are the relative injection 
efficiencies normalized to the baryon injection efficiency. The injection efficiencies tend to be higher 
for higher 70 shocks when the subshock is small. 
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shock speed if there is no subshock. Leptons with mass nie/mp = 0.01 may reach their maximum 
injection efficiency above a Lorentz factor of ^ 100. 
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addition, in the Fireball model large numbers of electron-positron pairs are produced. This suggests 
that the energy is carried predominately by the leptons due to their large particle number density. 
In this section, the characteristics of energy distributions are explored for various lepton to baryon 
particle number density ratios. 

Clearly, the energy carried by particles is in the ratio of their masses for particles with the 
same velocity. It was shown that when the leptons and baryons have equal particle densities, the 
baryons, and in particular the highly accelerated ones, carried significant energy and exerted the 
predominant backpressure on the shock to modify its velocity profile. As the lepton to baryon density 
ratio increases, the situation changes. Large numbers of leptons, compared to baryons, begin to exert 
their own backpressure on the shock velocity profile. Referring to Figure d. Ill for a lepton density 
a little higher than the baryon density, say for a case where n^jn^ ~ 10, the shock velocity profile 
is still dominated by baryon momentum and the smoothing takes place on the scale of the baryon's 
gyroradius. However, when leptons dominate, as is the case of a lepton to baryon density ratio of 
10"* and a lepton mass of O.Olmp, the smoothing results more from the backpressure of the energetic 
leptons and on the scale of the lepton gyroradius, recalling that smoothing always occurs in such a 
way as to conserve momentum and energy flux across the shock. The resulting smoothed shock has 
a shock width very small on the scale of the baryon diffusion length, and is on the order of the lepton 
diffusion length. Leptons now interact with the shock as baryons did when baryon momentum was 
dominant. Clearly, all baryons returning from downstream see the entire shock width. 

In Figure 16.81 the relative injection efficiency was defined at the onset of the power law for 
baryons, and a critical momentum po was used to define the lepton injection efficiency because the 
lepton power law started at approximately that point. In other words, few leptons received enough 
initial momentum after the first crossing to return for additional acceleration and develop the power 
law. The leptons that remained had to undergo repeated accelerations to develop enough momentum 
with their related diffusion length to see the entire shock, hence the momentum flux of leptons in the 
power law region was very low, due to both their small mass and relatively few particles remaining, 
compared to baryons. This definition of (relative) injection efficiency is somewhat misleading when 
the shock front becomes narrow on the scale length of most particles, and this situation occurs when 
the lepton to baryon number density ratio increases beyond 10^. As leptons begin to have more 
influence over the shock modification and carry more total energy, and as the shock sharpens as 
shown in Figure IH. Ill far more leptons can enter into the shock acceleration process after the first 
shock crossing. In other words, after the first shock crossing, the diffusion length of a much larger 
population of leptons extends over the entire lepton-modified shock width and the leptons start their 
power law, as shown in Figure 16.121 at point a. Point a is determined by the ratio of the masses 
because after the first crossing, as discussed earlier, the average speed for all particles is Urci and the 
momentum of a particle is 7rci'7iWrci. 

All baryons have a sufficient diffusion length to see the entire shock and a large population of 
leptons can now see the entire shock. Therefore, both baryons and leptons have the same power law 
dependence after the first crossing and the relation, Pa/pb = We/jTip results, as shown at the top of 
Figure 16.121 This point becomes well defined and takes the shape of the baryon injection efficiency 
plot as the lepton to baryon density ratio continues to increase. The initial p&rticle injection efficiency 
(i.e., in those cases where the power law begins after the first shock crossing) is defined as points 
a and 6. Defining the initial injection efficiency in this way indicates that the lepton and baryon 
particle injection efficiencies are nearly equal in Figure 16.121 but they occur at different momenta 
as mentioned above. In the previous section the mass difference was tied to the maximum relative 
injection efficiency through equation (|6.2|l . Now the initial injection efficiency is connected to the 
maximum relative injection efficiency through this equation. 

In order to make comparisons between lepton and baryon flux characteristics, they need to be 
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Figure 6.11: Velocity profiles steepen with increased lepton/baryon particle density ratios. The 
profile for a lepton density of 10* cm~^ with a lepton mass of O.Olmp is also smoothed, but on the 
scale of the lepton gyroradius and therefore it appears very sharp to baryons, whose gyroradius is 
used here for units on the horizontal axis. 
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Figure 6.12: Initial injection efficiencies are shown at a for leptons and at b for baryons after the 
first shock crossing of each species in a scenario where leptons dominate baryons by 10^. Due to 
the narrow or sharp shock, modified primarily by leptons, the injection efficiencies of all particles 
is relatively high and leptons achieve a power law after the first crossing, similar to baryons. Note 
that the baryon data were normalized by multiplying by the Ue/up ratio. 
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measured at the same momentum, therefore the relative injection efficiency wih be used as before, 
by taking the differences in logarithms between the lepton and baryon curves along the critical 
momentum line at pq (i.e., the momentum at point b in Figure l!).12|l . This was done for a range 
of lepton to baryon particle number densities and for a range of shock speeds, using a lepton to 
baryon mass ratio of 0.01. The results are expressed in Figure 15. 131 as the ratio of lepton to baryon 
energy flux as a function of lepton to baryon particle number density. When the particle number 
density is low, flux characteristics are dominated by baryons and the baryon-smoothed shock. As 
the leptons become more important and the shock smoothing begins to be influenced by leptons as 
well as baryons, leptons begin to share more of the energy density, and eventually there will be a 
point where the leptons and baryons share energy equally. This point, on the equipartition line in 
Figure 15.131 appears to be at approximately n^^jn^ 1600 for a lepton mass of O.OlWp. In other 
words, as the lepton to baryon particle number density increases from ~ 10^ to ~ 10'*, the spectral 
characteristics of shock accelerated baryons and leptons with mass — O.Olmp, and the shock 
profile itself, makes a transition from those dominated by baryons to those dominated by leptons. 
At higher number density ratios, both leptons and baryons behave the same. Leptons dominate 
completely above a number density ratio of about 10"* and, from that point on, the energy flux is 
simply proportional to the number density, hence the straight line (with a slope of 1). 

Additional energy fluxes were found at a number density ratio of 10^ for lepton masses of O.lmp 
and O.OOlmp, and from these a projection was possible down to the electron mass of (l/1836)TOp. 
Straight lines, each with a slope of 1, can be projected back to the equipartition line and, if this is 
done for the electron mass (designated as the solid star in Figure 15. 131 it can be estimated that the 
equipartition point for particles having the rest mass of an electron is at a lepton to baryon particle 
number density ratio of 3.2 x 10^. The measurements were collected from shocks with Lorentz factors 
of 2 and 12, the range of shock speeds in this study, and the results appear to be independent of 
these shock speeds above the equipartition line. 

These concepts are expressed in a different way in Figure 16.141 Here, fractional total energy 
above momentum p is shown for leptons and baryons for the case of relatively small lepton to baryon 
density in the top frame, and for a large lepton to baryon density in the lower frame. In the top 
frame, baryons carry the bulk of the total energy, even though lower in particle number because they 
entered into the acceleration process early and have a much higher density of accelerated particles as 
a result of baryons still dominating the shock smoothing. In the lower frame, leptons dominate the 
shock smoothing, also enter the acceleration process early, have a much higher density of accelerated 
particles (compared to leptons in the top frame) and, by their sheer numbers, now carry the bulk of 
the energy. In Figure 15.141 the curves are not normalized because the point here is to emphasized 
the actual fraction of total energy that is associated with each particle type. 

Noting the vertical dotted line in Figure 16.141 another parameter can be defined as energy 
efficiency (sometimes called acceleration efficiency as it was defined in section 5.3.3), similar to what 
was done in Figure l5.1()l Energy efficiency is the fraction of total energy remaining above momentum 
p. Since the dotted line in Figure 15. 141 is chosen, as before, at the onset of the power law portion 
of the baryon curve, this momentum is also the critical momentum, pa. Hence, energy efficiency (or 
acceleration efficiency) is 



Using this definition, ratios of the energy efficiencies can be measured for various lepton baryon 
density ratios and for various shock speeds, similar to what was done for injection efficiencies. These 
results were compiled for a lepton mass of 0.01. When they arc plotted on the same graph in Figure 
16.131 they overlay the me ~ 0.01 curve exactly, and this is not surprising. Injection efficiency and 



eacc(p > Po) 



fraction of total kinetic energy density above pq 
total kinetic energy density 



(6.3) 
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m„ = 0.01 line 
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Log.Q[lepton/baryon Number density] 



Figure 6.13: Partitioning of energy among leptons and baryons as a function of the relative number 

density of leptons. The energy density is split equally between leptons of mass O.Olmp and baryons 
when lepton density exceeds baryon density by approximately 1600. Leptons with the rest mass of 
the electron reach equipartition at a density ratio of 3.2 x 10^. 
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Figure 6.14: Fraction of total kinetic energy carried by baryons and leptons is shown for two 
lepton to baryon density ratios. In the top frame where the lepton density is 100, the kinetic energy 
carried by Icptorrs is up to four orders of magnitude below the faction carried by baryons. When 
the lepton density increase to 10^, as shown in the lower frame, the fraction of energy carried by 
leptons dominates. 



Chapter 6. Characteristics of lepton acceleration in relativistic parallel shocks 



82 



energy efficiency are closely related because injection efficiency is concerned with the fraction of 
particles that get injected into the acceleration process and energy efficiency is concerned with the 
fraction of kinetic energy carried by these same particles. Since they were defined at the point where 
each curve cnitcrcd the (constant) power law dependence, the most they could vary is by a constant. 
The constant is 1. 

6.4 Summary 

Leptons with various masses and baryons were subjected to a range of plane parallel modified shocks 
having Lorentz factors from 2 up to 20. The energy and momentum fluxes of these particles were 
studied as a function of unmodified subshock size and particle number density to try to better 
understand the parameters required for leptons to carry a significant portion of the energy in a 
gamma-ray burst. 

Assuming equal numbers of leptons and baryons, baryons dominate the shock smoothing pro- 
cess and lepton injection efficiencies are very sensitive to the size of any embedded subshock, the 
sensitivity increasing with decreasing lepton mass. 

For subshock fractional sizes, fs, above 0.4, the injection efficiencies of leptons appear to be 
insensitive to shock speeds over the range used in this study. In the case of negligible subshock, 
the injection efficiencies of leptons are sensitive to the shock speed, and projections suggest that for 
shock Lorentz factors on the order of 100, the shock front may be sufficiently narrow to increase 
the relative injection efficiency of leptons with mass O.Olrrip to a maximum of 10~^'^ of the baryon 
injection efficiency. A lepton mass of nie ^ (l/1836)mp may require a shock Lorentz factor greater 
than 1000 to reach a maximum relative injection efficiency of about 10~^'^ that of baryons. The 
overall conclusion is that, given the assumptions in the Monte Carlo model, leptons can never carry 
enough energy to explain the spectra observed in gamma-ray bursts if the lepton to baryon particle 
densities are of the same order, even for ultrarelativistic shock speeds. 

The injection efficiencies of leptons and baryons were studied in the case where the lepton to 
baryon number density ratios were varied and a lepton to baryon mass ratio of 0.01 was used. It 
was found that, as expected, for low number density ratios, baryons dominate the shock smoothing 
and the lepton injection efficiencies and energy efficiencies remained low. As the number density 
ratio increased to the 10^ range, leptons began to have a strong influence on the modification of the 
shock velocity profile, the shock became much sharper, and a much larger population of leptons were 
injected into the shock process. The energy equipartition point of leptons with mass of O.OlTOp was 
at the particle number density ratio of 1600. Above a number density ratio of 10^, leptons were found 
to dominate the shock smoothing process and the energy flux was directly proportional to number 
density ratio. Using this method, an extrapolation was made down to the mass of the electron, and 
it was determined that the energy equipartition point for electrons (and positrons) in plane parallel 
shocks is at a lepton to baryon number density ratio of 3.2 x 10"''. This result was independent of 
shock speeds with Lorentz factors of 2 to 12 used in this study. However, the energy equipartition 
point may be independent of shocks with any higher Lorentz factor for the following reason. When 
the lepton to baryon partic;le densities are of the same order, the maximum lepton relative injection 
efficiency occurs at ultrarelativistic shock speeds and it is far below baryon injection efficiency (at 
the same momentum). Hence, the only way the lepton efficiency can dominate, under the same 
acceleration assumptions, is by increasing the lepton to baryon number density ratio. When the 
lepton density ratio is high enough that leptons dominate the shock modification process, shock 
speed no longer plays a roll in increasing injection efficiency. The high number density ratio of 
leptons to baryons is required for two reasons. The low mass of the leptons means that each lepton 
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carries far less energy than baryons, in proportion to their masses at the same speeds, therefore large 
numbers of leptons are necessary to carry an equal portion of the total energy. More importantly, 
there must be a sufficient lepton density to dominate the shock smoothing so that larger numbers 
of leptons can be injected into the shock acceleration process. 

The primary finding in this study was that for plane modified shocks with the magnetic field 
parallel to the shock normal, and with shock Lorentz factors in the 2 to 12 range, the particle number 
density ratio between leptons and baryons must be above 3.2 x 10^ for leptons with the rest mass 
of the electron to carry sufficient kinetic energy in gamma-ray bursts and account for the observed 
GRB spectra. The energy equipartition point appears to be independent of shock speed. 



Chapter 7 



Conclusions 



Chapter summaries can be found at the end of each chapter, but the accomphshments and main 
conclusions of the entire dissertation are summarized here. Despite the fact that relativistic shock 
theory has concentrated almost exclusively on test-particle acceleration, it is likely that relativistic 
shocks are not test-particle but inject and accelerate particles efRciently. The reason is that regardless 
of the ambient far upstream conditions, all particles that are overtaken by a relativistic shock will 
receive a large boost in energy /S.E ^ 7ioimc^ in their first shock crossing. Thus, virtually all of 
the particles in the downstream region of an unmodified shock are strongly relativistic with u ~ c. 
The ability to overtake the shock from downstream and be further accelerated depends only on 
the particle speed (Eq. 15. 7|) , the presence of magnetic waves or turbulence with sufficient power in 
wavelengths on the order of the particle gyroradii to isotropize the downstream distributions, and the 
pitch angle of the returning particle fEa. I5.4|l . It is generally assumed that the necessary magnetic 
turbulence is self-generated and if enough turbulence is generated to scatter high momentum particles 
(with very low densities) that constitute a test-particle power law, there should be enough generated 
to isotropize lower momentum particles (which carry the bulk of the density). If acceleration can 
occur at all, it is likely to occur efRciently making it necessary to calculate the shock structure and 
particle acceleration self-consistently. Furthermore, if relativistic shock theory is to be applied to 
gamma-ray bursts, where high conversion efficiencies are generally assumed, nonlinear effects must 
be calculated. 

When energetic particles are generated in sufficient numbers, the conservation of momentum 
and energy requires that their backpressure modify the shock structure. Two basic effects occur: 
a precursor is formed when the upstream plasma is slowed by the backpressure of the accelerated 
particles and the overall compression ratio changes from the test-particle value as a result of high 
energy particles escaping and/or a change in the shocked plasma's adiabatic index T. As indicated 
by the 70 = 1.4 example (Section 5.3.1), mildly relativistic shocks act as nonrelativistic ones showing 
a dramatic weakening of the subshock combined with a large increase in r ('Fig. l5.injl . These changes 
result in a particle distribution which is both steeper than the test-particle power law at low momenta 
and flatter at high momenta fFig. 15.11)) . 

In faster shocks (i.e., 70 ^ 3), the initial test-particle spectrum is steep enough that particle 
escape is unimportant so only changes in F determine r (Ea. I5.l5|l . In contrast to nonrelativistic 
shocks where the production of relativistic particles causes the compression ratio to increase, it is 
shown that r decreases smoothly to 3 as 70 increases and the fraction of fully relativistic shocked 
particles approaches one (Fig. I5.H|I . 

An important result is that efficient, mildly relativistic shocks do not produce particle spectra 
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close to the so-called 'universal' power law having a ^ 4.23. This may be important for the inter- 
pretation of gamma-ray bursts since the internal shocks assumed responsible for converting the bulk 
kinetic energy of the fireball into internal particle energy may be mildly relativistic and the external 
shocks, believed responsible for producin g gamma-ray burst afterglows, will always go through a 
mildly relativistic phase fsee IPiranl Il999l for a comprehensive review of gamma-ray bursts). 

7.1 Relativistic magnetohydrodynamic Monte Carlo model 

The basic shock acceleration physics is independent of the speed of the shock, but the mathematical 
modeling of the shock process depends on whether or not particle speeds are large compared to the 
shock speed. The relativistic Monte Carlo model used for this study, allows particle speeds to be 
comparable to shock speeds and, aside from computational limits, the model allows shocks to have 
arbitrary Lorentz factors. 

The fully relativistic nonlinear Monte Carlo model, under development from a nonrelativistic 
model, will simulate kinematic particle acceleration in steady-state modified shocks with a gyrotropic 
pressure tensor and magnetic fields at oblique angles when it is completed. In its present state, 
operating with the magnetic field parallel to the shock normal, the shocks can be self-consistently 
modified while conserving particles, momentum, and energy across the shock, and while allowing 
charged particles to gain energy by first-order Fermi acceleration, and escape when they reach a 
high enough energy. The model can also be used with oblique angles and anisotropic pressure for 
unmodified test particle shocks. 

New flux conservation laws of energy and momentum were developed by setting the divergence 
of the stress-energy tensor equal to zero and applying the 4-dimensional version of Gauss' theorem. 
Combining these equations with Maxwell equations yields a completely general set of relativistic 
magnetohydrodynamic jump conditions across the shock. A new equation of state was developed 
from the conservation of total energy that incorporates oblique magnetic fields and gyrotropic pres- 
sure. A new equation was developed from kinetic theory that gives an excellent approximation to 
the adiabatic index at intermediate shock speeds. The adiabatic index becomes the normal ratio of 
specific heats at the nonrelativistic limit (5/3) and at the ultrarelativistic limit (4/3). An anisotropy 
parameter was introduced for the gyrotropic pressure tensor to allow an analytic solution to the equa- 
tions. The Monte Carlo code was further developed with new relativistic momentum transformation 
equations and a method of handling fluxes at grid zones for relativistic shock modification. 

The Monte Carlo model was used first to duplicate well known test particle power law results 
in fully nonrelativistic and ultrarelativistic parallel shocks and was shown to move all parameters 
smoothly between these two shock speed regimes. The model was used in the trans-relativistic regime 
and can self-consistently determine compression ratios and adiabatic indices in close agreement 
with recent analytic results of Kirk et al. (2000). A number of solutions were presented where the 
downstream pressure is anisotropic. This is believed to be the first presentation of oblique shock 
solutions which apply smoothly for nonrelativistic, trans-relativistic, and ultrarelativistic shocks. 

The results presented in Chapter 5 assumed that no first-order Fermi acceleration occurs, but 
they apply directly to test-particle acceleration where the energy density in accelerated particles is 
small. In that case, the differences in the compression ratio, r, produced by large magnetic fields 
and/or anisotropic pressures, even at ultrarelativistic speeds, become important since changes in 
r map directly to changes in the power law index of the accelerated spectrum. This power law is 
the most important characteristic of test-particle Fermi acceleration and is often associated with 
ultrarelativistic internal shocks in gamma-ray burst (GRB) models. In standard GRB models, the 
rapidly expanding fireball cools, converting the internal energy of the hot plasma into kinetic energy 
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of the relativistically moving ejecta and electron-positron pairs. At the point where the fireball 
becomes optically thin and the GRB we see is emitted, the matter is too cool to emit gamma- 
rays unless some mechanism can efficiently re-convert the kinetic energy back into random internal 
energy, i.e., unless some particle acc eleration process takes p l ace. Fermi sh ock acceleration is widely 
believed to be this mechanism fe.g..'Rees fc Meszarosl fl992l: IPira 

Internal flows in GRBs may be ultrarelativistic, but when a fast flow overtakes a slower flow, the 
resulting shocks may be in the trans-relativistic range and may be responsible for the early transient 
gamma-ray intensity peaks. The outer blast wave shock slows as it sweeps up and accelerates 
interstellar material. This shock, believed to produce long-lasting afterglows, will eventually always 
go through a trans-relativistic phase, and the results apply to particle acceleration here as well. 



7.2 Application to lepton and baryon acceleration 

The Monte Carlo model was applied to lepton and baryon particle acceleration for the case of 
moderately relativistic, smoothed shocks with the shock normal parallel to the magnetic field. Using 
leptons of various light masses, and equal particle densities, it was found that unmodified subshocks 
of various fractional sizes within the main modified shock produced a large change in the lepton 
flux. The sensitivity to the larger fractional sizes of the subshock was little affected by the speed of 
the shock in the range of Lorentz factors between 2 and 12. However, in the case of no subshock, or 
small subshock, it was found that the lepton to baryon relative injection efficiency ratio, based on 
the critical momentum, po, where baryons start the power law distribution, increased substantially 
as the shock speed increased. For leptons with mass of O.OImp that of the baryon, leptons reach 
their maximum relative injection efficiency of approximately 10^^'^ that of baryons around a shock 
Lorentz factor of 100. For a lepton mass of rrie ^ (I/I836)mp, the shock Lorentz factor may need to 
be greater than 1000 to reach a maximum relative injection efhciency of about lO^'*-^ that of baryons. 
Such speeds may briefly occur in the initial outflow of a G RB, but it w ould require the entrained total 
baryonic mass to be very smaU (i.e., less than 10-*Mq; Esii (200^. The conclusion here is that, 
given the assumptions of this model with particles energized by diffusive shock acceleration, leptons 
cannot carry sufficient energy to explain the observed spectra if the lepton and baryon densities are 
of the same order. 

Present gamma-ray burst models require leptons to carry a significant fraction of the total energy 
density of the burst. The most likely possibility is that there are far more leptons than baryons. 
Arguments connected with opacity and variability timescales, etc fe.g.. iPira require highly 

relativistic outflow in GRBs and this is only possible if the total baryonic mass is low. If the total 
baryonic mass is low and if electron-positron pairs are created in large numbers by the init i al gam ma- 
ray photon burst while the burst is still opaque ("r a mirez- Ruiz. MacFadven. fc LazzatA l2002l) and 
by neutrinos if the baryonic density is low (Rees, 2000), then lepton to baryon particle density ratio 
should be high. If it is high enough, then leptons can carry the bulk of the available burst energy. 

The question then is what ratio of leptons to baryons would be required for equal sharing of the 
energy between the two species. The characteristics of leptons and baryons with various particle 
number densities were studied for plane parallel, relativistic, modified shocks in the same speed 
range. Initial energy and momentum of the particles, after the first shock crossing, was found to 
be in the ratio of their masses as expected. For moderately high lepton to baryon number density 
ratios (e.g., 10 - 100), the smoothed shock velocity profile modification was completely dominated 
by baryons. Using a lepton mass of O.OImp, it was found that the shock and spectral characteristics 
make a transition from baryon dominated to lepton dominated in the lepton to baryon number 
density ratio rang e of 10^ to 10"*. Leptons with mass O.OlTOp shared energy density equally with 
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baryons when the particle number density ratio was approximately 1600. Above this point, leptons 
begin to completely dominate the energy and energy flux, and was found to be directly proportional 
to number density ratio. 

By using leptons of various masses, a projection was possible down to the mass of the electron 
{me/mp ^ 1/1836). An important and new result was that the energy density equipartition point 
between electron-positron pairs and baryons was at a lepton to baryon particle number density 
ratio of 3.2 x 10^, and it appears to be independent of the shock Lorentz factor. This suggests 
the possibility that, over a wide range of shock conditions, the number density ratio of leptons to 
baryons required for equipartition of kinetic energy in a GRB may stay roughly constant, similar to 
the constant value found above for parallel shocks. 

After the Monte Carlo model is completed, the model will be used to further investigate lepton 
and baryon acceleration using oblique nonlinear shocks. This may better establish the lepton to 
baryon number density ratio required for equipartition. 
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Definitions and Fundamental Relationships 



Alfven velocity va '^here B is tlie magnetic field strength, 7 is the pertinent Lorentz 



factor and p is the rest mass density. The Alfven velocity refers to the speed of propagation 
of a transverse magnetic disturbance along the direction of the magnetic field. 

Gyroradius — for charged particles in the presence of a magnetic field, where p is the 

particle's momentum, q is the unit electric charge, Zq is the particle's total charge {Z is the 
total charge number), and B is the magnetic field strength. 

Mean free path A — rjrg, is the average distance between collisions, if the particles undergo colli- 
sions. For coUisionZess interactions (e.g., through magnetic field tubulence), a similar concept 
can be defined as the average distance required for the particle to change direction by 90 
degrees. Bohm diffusion is assumed when the gyrofactor r] — I. Higher values of rj imply 
smoother magnetic fields, less scattering, and longer mean free paths. 

Diffusion coefficient n = •^Xv = irjroi) = i^^rr = ft— , where 3 is the speed of the 

3 3 3 ZqB 3 ZqB 

particle. 

Diffusion length L = — = 77 — ~ 77 ^ (x pv, where un is the shock speed. 
Wo 3 Wo 3 uo 

Collision time '■c = ~ = is a means of estimating the time for a particle to scatter through 
an angle of 90 degrees. 



Debye length Ad cx y ~^ ^ measure of the electrostatic shielding distance around a charged 

particle in a plasma l|Shi]L Il993) . E refers to the kinetic energy of electrons, n is the total 
particle number density, and q is the unit electric charge. 

... "DC 

Rigidity R ~ provides a measure of the sensitivity of the momentum of a particle to a magnetic 
Zq 

field. The higher the rigidity, the less influence a magnetic field has on the path of the particle. 

Energy and momentum flux distributions distributions f{p) or g{E) are defined here as the 
total number of passages of particles through a grid plane from any direction. Referring to 
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Figure A.l: Momentum space diagram showing an element of momentum space (Pp 



Figure IA.1L the total number of particles in momentum space is 



f{p)d-'p = At: / p'f{p)dp 
1 Jo 

provided the integration takes place where the distribution is isotropic; then 



Airp^ f{p)dp gives the 



# of particles 

z::^ 



in dp 



In the same way, if 



(A.l) 



(A.2) 



(A.3) 



dJ 



represents the total number of particles in energy space, where represents the differential 
flux in particles/ (cm^ • second • ster • E), then 



# of particles 



cm 



in dE 



47r dJ ,„ . 
— -jT^dbj gives the 
V dh ° 

and the densities can be equated, yielding: 

dE '^P-'^P'dE 

Finally, using the relation E'^ ~ p^c? + {mc^Y , equation (|A.5p simplifies to 



(A.4) 



(A.5) 



(A.6) 
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Determining the ratio of specific heats 



B.l Kinetic pressure and the equation of state 

The total distribution of particle momenta in a plasma gives the plasma a certain internal energy 
which has associated with it an adiabatic exponent or index. The discussion that follows develops a 
relationship between this adiabatic index and the calculated average momentum, or kinetic energy 
of the particles, based on simple statistical physics. The relationship can be used to estimate the 
adiabatic index in the downstream plasma frame of a relativistic shock. The concepts are well known 
and versions of the derivation belo w are found in most basic texts that discuss kinetic theory (e.g., 
iFevnman. Leighton fc SandslllQS'^ . This overview is for the convenience of the reader. 

The adiabatic equation of state for an ideal gas of point particles, relating the pressure and 
volume of a closed system, may be written as 

PV^ = C (B.l) 

where P is the kinetic pressure. Vis the volume under consideration, F is the adiabatic index, and C is 
a proportionality constant. The significance of point particles is that the particles have no rotational 
or vibrational contributions and all of the kinetic energy of the particles is due to translation. 
Taking logarithms and then the derivatives in the equation above, one has 

^ + F^ = or VdP + TPdV = (B.2) 



Adding PdV to both sides and rearranging, 

PdV + VdP 
F- 1 



= -PdV = dU. (B.3) 



Since —PdV represents an adiabatic incremental compression of a volume at constant pressure 
and dU is the corresponding increase in total internal energy, and since the left hand side has the 
total derivative of PV, one can integrate, with the appropriate limits, and write: 

PV UP , , 

1—1=^ or y-u^^^. (B.4) 

where u is the total internal energy density of the particle ensemble. 

The exponent, F, is sometimes called the ratio of specific heats, but this is only accurate in 
the nonrelativistic (F = 5/3) and ultrarelativistic (F — 4/3) limits. Between these two limits, F is 
simply the exponent for the adiabatic relation or the adiabatic index. 
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The total energy density e and isotropic pressure P are related through a combination of the 
adiabatic equation of state (i.e., equation ljB.4|) and the conservation of energy, which forms another 
equation of state: 

e = + pc^ (B.5) 

B.2 Nonrelativistic case 

Consider a rigid rectangular chamber with a movable wall on the right end as shown in Figure 
IB. II The chamber contains ideal point particles of particle density n. The volume of the chamber 
has dimensions x, y, z and the movable wall with area A — yz on the right end is pushed into 
the volume a distance —dx by a normal force F. Then the reversible work done on the volume is 
dW — F{—dx) — PA{—dx) — —PdV, where —dV is the corresponding decrease in volume and P is 
the pressure inside the box. 




Figure B.l: A rectangular volume containing ideal particles under adiabatic compression from a 
force F. 

The total momentum change by a particle moving in the x direction inside the volume and 
elastically colliding with the movable wall is 5p = 2px. The total number of particles that collide 
with the movable wall in time 5t is nVxStA. Hence, the force on the movable wall inside the volume is 

F = = {nVxA){2px) and the pressure is P = 2npxVx- The particle velocities have a wide range 

of speeds and directions and, in particular, only half the particles in the x direction move toward 
the movable wall. Taking these considerations into account while averaging over the velocities, the 
pressure becomes P = n < p^v^ >■ If there are no magnetic fields and the x,y and z velocity 

distributions are equivalent, then < PxVx >= ^ < p- v >, and the isotropic pressure can be written 
as 

P = < p - V > (B.6) 
It might be noted here, that if magnetic fields were present, one might take averages of the momenta 
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parallel and perpendicular to the magnetic field, assuming rotational symmetry about the field vec- 
tor. In that case one could have a gyrotropic pressure (tensor) with the following space components 
along the diagonal: 

P» = n <p - V >|| (B.7) 



and 



In general, the momentum for a particle with mass m and speed v is 

JeI^ - {mc^ f 

p = — = mc\/ 7^ — 1 = •mc\/ ^"^(P = 77711; (B-9) 

where the Lorentz factor 7 = (1 j) 2 . 

Using this relationship for momentum one can write, for isotropic pressure, 

P = '^n{-fmv^) . (B.IO) 
For speeds where 7 ~ 1, equation (|2.6(l can be approximated to first order as 

where U is the total internal energy and is the total kinetic energy of ideal point particles in the 
volume under consideration. 

2 



Comparing equation ljB.ll() to equation ljB.4(l above, it can be seen that F — 1 corresponds to ^ , 
5 

which yields F = ^ . This is the value for the adiabatic ratio of specific heats for ideal point particles 

having nonrelativistic speeds where 7 ~ 1. 

An error is introduced if the expression ^mv^ /2 is interpreted to be kinetic energy at mildly 
relativistic speeds where the appropriate expression for kinetic energy should be 7710^(7 — 1). For 
example, a 0.1 per cent error occurs, and will continue to increase for v ^ 0.1c. At 1; ~ 0.6c, the 
error increases to about five per cent. 

B.3 Ultrarelativistic case 

The same ideas can be used to derive an expression for isotropic pressure for ultrarelativistic particles 
in the volume described above, and the resulting relationship will be the same as in equation ljB.6|l . 



i.e. P — i?! < p - V >. However, for true ultrarelativistic particles, the momentum of the particles is 
p = — kinetic energy Ek ~ Etot — "yrnc^ and v ^ c. The equation for isotropic pressure becomes: 

P ^ \n (jmc^) ^ \^ (B.12) 



1 , 4 



and by the same comparison, this time equation (|B.12p to equation l|B.4p . F — 1 = ^ ^^'^ "'^ ~ 3' 
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An error is introduced if the ultrarelativistic case is assumed at less relativistic speeds where 
the rest mass cannot quite be neglected. For example, a 0.1 per cent error occurs, and will continue 
to slowly increase for Lorentz factors 7 ^ 500, but at 7 ~ 10, the error has only increased to 
approximately five per cent. 



Bibliography 



Achterberg, A., Gallant, Y. A., Kirk, J. G., & Guthmann, A. W., Particle acceleration by ultra- 
relativistic shocks: theory and simulations, 2001, M.N.R.A.S., 328, 393 

Appl, S. & Camenzind, M., Shock conditions for relativistic MHD jets, 1988, A&A, 206, 258 

Axford, W. I., Leer, E., & Skadron, G., The Acceleration of Cosmic Rays by Shock Waves, 1977, in 
Proc. 15th ICRC (Plovdiv), 11, 132 

Ballard, K.R. & Heavens, A.F., First-order Fermi Acceleration at Oblique Relativistic Magnetohy- 
dromagnetic Shocks, 1991, M.N.R.A.S., 251, 438 

Bednarz, J., & Ostrowski, M., The acceleration time-scale for first- order Fermi acceleration in rela- 
tivistic shock waves, 1996, M.N.R.A.S., 283, 447 

Bednarz, J. & Ostrowski, M., Energy Spectra of Cosmic Rays Accelerated at Ultrarelativistic Shock 
Waves, 1998, Phys. Rev. Letts, 80, 3911 

Bell, A. R., The Acceleration of Cosmic Rays in Shock Fronts - 1, 1978a, M.N.R.A.S., 182, 147-156 

Bell, A. R., The Acceleration of Cosmic Rays in Shock Fronts - II, 1978b, M.N.R.A.S., 182, 443-445 

Bell, A.R., The non-linear self-regulation of cosmic ray acceleration at shocks, 1987, M.N.R.A.S., 
225, 615-626 

Berczhko, E. G., & Ellison, D. C, A Simple Model of Nonlinear Diffusive Shock Acceleration, 1999, 
ApJ, 526, 385 

Berezhko, E. G., & Volk, H. J., Kinetic Theory of Cosmic Rays and Gamma Rays in Supernova 
Remnants. I. Uniform Interstellar Medium, 1997, Astroparticle Phys., 7, 183 

Blandford, R. D., & Eichler, D, Cosmic Ray Acceleration in Astrophysics, 1987, Phys. Repts., 154, 
1-75 

Blandford, R. D. & McKee, C. F., Fluid Dynamics of Relativistic Blast Waves, 1976, Phys. Fluids, 
19, No.8, 1130 

Blandford, R.D., & Ostriker, J. P., Particle Acceleration by Astrophysical Shocks, 1978, ApJ, 221, 
L29 

Cabannes, H., Theoretical Magnetofluid Dynamics Academic Press (1970) 

Cheng, K.S. & Lu, T., Gamma-Ray Bursts: Afterglows and Central Engines, 2001, Chin. J. Astron. 
Astrophys., 1, No. 1, 1-20 



94 



BIBLIOGRAPHY 



95 



de Hoffman, F. & Teller, E., Magneto- Hydrodynamic Shocks, 1950, Phys. Rev., 80, 692 

Dieckmann, M.E., et.al., Electron Acceleration due to High Frequency Instabilities at Subpernova 
Remnant Shocks, 2000, A&A, 356, 377-388 

Drury, L. O'C, An Introduction to the Theory of Diffusive Shock Acceleration of Energetic Particles 
in Tenuous Plasmas, 1983, Rep. Prog. Phys., 46, 973-1027 

Dwarkadas, V.V. & Chevalier, R.A., Interaction of Type la Supernovae with their Surroundings, 
1998, ApJ, 497, 807-823 

Ellison, D.C., Monte Carlo Simulation of Collisionless Shock Acceleration, PhD. Thesis, 1981, 
Catholic University 

Ellison, D.C., Shock Acceleration of Diffuse Ions at the Earth's Bow Shock: Acceleration Efficiency 
andA/Z Enhancement, 1985, J.G.R., 90-Al, 29-38 

Ellison, D. C. 1991a, FermA Particle Acceleration in Relativistic Shocks: Spectra and Efficiencies from 
Modified Shocks, Proc. ICRR International Symposium on "Astrophysical Aspects of the Most 
Energetic Cosmic Rays," p. 281, Eds., M. Nagano and F. Takahara, Wold Scientific, Bangalore. 

Ellison, D. C. 1991b, Fermi Particle acceleration in relativistic shocks: Preliminary nonlinear results, 
in "Relativistic Hadrons in Cosmic Compact Objects," p. 101, Eds., A. A. Zdziarski and M. Sikora, 
Springer- Verlag, Berlin. 

Ellison, D. C, Baring, M. G., & Jones, F. C, Acceleration Rates and Injection Efficiencies in Oblique 

Shocks, 1995, ApJ, 453, 873-882 

Ellison, D. C, Baring, M. G., & Jones, F. C, Nonlinear Particle Acceleration in Oblique Shocks, 
1996, ApJ, 473, 1029 

Ellison, D.D. & Double, G.P., Nonlinear Particle Acceleration in Relativistic Shocks, 2002, As- 
troparticle Phys., 18-3, 213-228 

Ellison, D. C, & Eichler, D. Monte Carlo Shock-Like Solutions to the Boltzmann Equation With 
Collective Scattering, 1984, ApJ, 286, 691 

Ellison, D. C, Giacalonc, J., Burgess, D., and Schwartz, S. J., Simulations of Particle Acceleration 
in Parallel Shocks: Direct Comparison Between Monte Carlo and One- Dimensional Hybrid Codes, 
1993, J.G.R., 98, 21,085 

Ellison, D. C, Jones, F. C, & Baring, M. G., Direct Acceleration of Pickup Ions at the Solar Wind 
Termination Shock: The Production of Anomalous Cosmic Rays, 1999, ApJ, 512, 403 

Ellison, D.C., Jones, F.C., and Reynolds, S.P., First-Order Fermi Particle Acceleration at Relativistic 
Shocks, 1990, ApJ, 360, 702 

Ellison, D.C., Mobius, E., & Paschmann, G., Particle Injection and Acceleration at the Earth's Bow 
Shock: Comparison of Upstream, and Downstream Events, 1990, ApJ, 352, 376 

Ellison, D. C, & Reynolds, S. P., ^ determination of relativistic shock jump conditions using Monte 
Carlo techniques, 1991, ApJ, 378, 214 

Fermi, E., 1949, Phys. Rev., 75, 1169-74 (Collected Papers, vol II, pp.655-55) 



BIBLIOGRAPHY 



96 



Feynman, R.P., Leighton, R.B. & Sands, M., Feynman Lectures on Physics, 1963, Addison- Wesley 

Gallant, Y. A., Particle Acceleration at Relativistic Shocks, preprint, 2002, |astro-ph/0 201243| 

Gallant, Y. A. & Achterberg, A., Ultra-high- energy cosmic ray acceleration by relativistic blast waves, 
1999, M.N.R.A.S., 305, L6 

Gleeve, A. A., The Heating and Acceleration of Electrons by Shocks, 1984, Adv. Space Res., 4, 
No.2-3, 255-263 

Heavens, A. F. & Drury, L. O'C, Relativistic shocks and particle acceleration, 1988, M.N.R.A.S., 
235, 997 

Hollweg, J., On the Origin of Solar Spicules, 1982, ApJ, 257, 345-353 

Jackson, J.D., Classical Electrodynamics, 2nd ed., 1975, John Wiley & Sons, NY 

Jones, F.C. & Ellison, D.C., Noncoplanar Magnetic Fields, Shock Potentials, and Ion Deflection, 
1987, J.G.R., 92 No. AlO, 11205 

Jones, F.C., & Ellison, D.C., The Plasma Physics of Shock Acceleration, 1991, Space Sci. Rev., 58, 
259 

Kirk, J. G., 1988, Particle Acceleration at Relativistic Shock Fronts, Thesis, Dr. rer. nat. habil., 
Ludwig-Maximillians-Universitat, Miinchen, Germany. 

Kirk, J.G. & DufFy, P., Particle Acceleration and Relativistic Shocks, J. Phys. G: Nucl. Part. Phys., 
1999, 25, R163 

Kirk, J. G., Guthmann, A. W., Gallant, Y. A., Achterberg, A., Particle Acceleration at Ultrarela- 
tivistic Shocks: An Eigenfunction Method, 2000, ApJ, 542, 235 

Kirk, J. G. & Schneider, P., On the acceleration of charged particles at relativistic shock fronts, 
1987a, ApJ, 315, 425 

Kirk, J. G., & Schneider, P., Particle Acceleration at Shocks: A Monte Carlo Method, 1987b, ApJ, 
322, 256 

Kirk, J.G. & Webb, G.M., Cosmic-ray Hydrodynamics at Relativistic Shocks, 1988, ApJ, 331, 336 

Krymsky,G.F., 1977, Dokl. Akad. Nauk SSSR, 234, 1306 (Engl, transl. Sov. Phys.-Dokl., 23, 327), 
1981, Izv. Akad. Nauk SSSR Ser. Fiz., 45, 461 

Landau, L.D. & Lifshitz, E.M., Fluid Mechanics, 1959, Adison- Wesley Pub. 

Landau, L.D. & Lifshitz, E.M., The Classical Theory of Fields, 1962, Adison- Wesley Pub. 

Levinson, A., Electron Injection in Collisionless Shocks, 1992, ApJ, 401, 73-80 

Levinson, A., On the Injection of Electrons in Oblique Shocks, 1996, M.N.R.A.S., 278, 1018-1024 

Lichnerowicz, A., 1967, Relativistic Hydrodynamics and Magnetohydrodynamics, Benjamin 

Malkov, M., 1998, Phys. Rev. E, 58, 4911 



BIBLIOGRAPHY 



97 



McClements, K.G., et.al., Acceleration of cosmic ray electrons by ion-excited waves at quasi- 
perpendicular shocks, 1997, M.N.R.A.S., 291, 241-249 

Meszaros, P. & Rccs, M.J., Poynting Jets from Black Holes and Cosmological Gamma-ray Bursts, 
1997, ApJ, 482, L29-L32 

Michel, F.C., 1981, ApJ, 247, 664 

Ostrowski, M., Monte Carlo simulations of energetic particle transport in weakly inhomogeneous 
magnetic fields - I. Particle acceleration in relativistic shock waves with oblique magnetic fields, 
1991, M.N.R.A.S., 249, 551 

Peacock, J. A., Fermi acceleration by relativistic shock waves, 1981, M.N.R.A.S., 196, 135 

Pelletier, G., Cosmic ray acceleration and nonlinear relativistic wavefronts, 1999, A&A, 350, 705 

Pelletier, G., & Marcowith, A., Nonlinear dynamics in the relativistic plasma of astrophysical high- 
energy sources, 1998, ApJ, 502, 598 

Piran, T., 1997 Unsolved Problems in Astrophsics, ed. J.N. Bahcall & J. P. Ostriker, Princeton 
University Press 

Piran, T., Gamma-ray bursts and the fireball model, 1999, Phys. Repts., 314, 575 

Ramircz-Ruiz, E., MacFadycn, A. I., & Lazzati,D., Precursors and pair loading from erupting 
fireballs 2002, M.N.R.A.S., 331, 197-202 

Rees, M.J., A Review of Gamma Ray Bursts, 2000, Nuc. Phys. A, A663&:664, 42c-55c 

Rees, M.J., & , Meszaros, P. Relativistic fireballs - Energy conversion and time-scales 1992, 
M.N.R.A.S., 258, 41 

Rybicki, G.B. & Lightman, A. P., Radiative Processes in Astrophysics, 1979, John Wiley & Sons 
(ISBN 0-471-82759-2) 

Scheuer, P.A.G. (1987) Astrophysical Jets and Their Engines, ed. W. Kundt, NATO ASI SEries, 
Reidel Dordrecht, 208, 129 

Schlickeiser, R. & Dermer, C.D., Proton and electron acceleration through magnetic turbulence in 
relativistic outflows, 2000, A&A, 551, 946-972 

Schmitz, H & Chapman, S.C., Electron Preacceleration Mechanisms in the Foot Region of High 
Alfvemc Mach Number Shocks, 2002, ApJ, 579, 327-336 

Schneider, P. & Kirk, J. G., Fermi acceleration at shocks with arbitrary velocity profiles, 1987, ApJ, 
323, L87 

Shu, F.H., 1992, The Physics of Astrophysics, Vol. H: Gas Dynamics, University Science Books 

Tan, J.C., Matzner, C.D. & McKee, C.F., Trans-relativistic blast waves in supernovae as gamma-ray 
burst progenitors, 2001, ApJ, 551, 946-972 

Taub, A.H., Relativistic Rankine-Hugoniot Equations, 1948, Phys. Rev., 74, 328-334 

Tidman, D.A. & Krall, N.A. (1971), Shock waves in collisionless plasmas, Wiley Science 



BIBLIOGRAPHY 



98 



Tolman, R.C., Relativity, Thermodynamics and Cosmology, University Press, Oxford (1934); 
unabridged Dover paperback reprint (1987) 

Vietri, M., The acceleration of ultra-high-energy cosmic rays in gamma-ray bursts, 1995, ApJ, 453, 
883 

Webb, G.M., Zank, G.P., & MccKenzie, J.F., Relativistic Oblique Magnetohydrodynamic Shocks, 
1987, J. Plasma Rhys., 37, part 1, 117 

Zel'dovich, Y.B. & Raiser, Y.P., Physics of Shock Waves and High- Temperature Hydrodynamic 
Phenomena, Academic Press, 1967 



